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F po( e ) 

g 

ge 

gL 


dummy integration variable for area 

surface area between the mass of material compacted by the landing 
surface and the mass of uncrushed crushable material 

horizontal planar projection of A c i 

surface area of lower portion of payload over which shear and/or 
normal stresses act to decelerate payload 

horizontal planar projection of Ap 0 and, in the absence of shear 
effects, Api 

surface area between the mass of material compacted by the payload 
and the mass of uncrushed crushable material 

horizontal cross-sectional area in the volume V c of figures 3 
and 4 (the volume that would lie beneath the landing surface if 
there had been no crushing) 

dimen sional d ummy integration variable used in equation (A27) ; 


dimen sional dum my integration variable used in equation (A29) ; 

= yp 2 i ~ 2 

dimensionless variable describing the relative motion of payload 

9p “ 9 

penetration; defined as — — - — in equation (A39) 

e k R P 

maximum value of e during payload penetration 

Rp 

identical to F cl (z) if — - y is substituted for z 

dimensionless vertical crushing force acting over area A cl and 
tending to decelerate the mass above A cl ; defined in equation 
(A 36) 

dimensionless vertical crushing force acting over area Ap 1 and 
tending to decelerate the mass above A pl ; defined in equation 
(A45) 

general symbol for acceleration due to gravity 

value of g on earth (32.17 ft/sec 2 herein) 

value of g at the landing site (12.3 ft/sec 2 herein) 
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local height of the volume which would lie beneath the landing sur- 
face if there had been no crushing; the local height of V c ; see 
figures 2, 3, and 4 

local height of the compacted volume of material compacted by the 
landing surface; the local height of m cl ; defined in equation 
(A12) ; see figures 2, 3, and 4 


local height of the compacted volume of material compacted by the 
payload; the local height of mpi; defined in equation (A12) ; see 
figures 2, 3, and 4 


defined in equation (17) as 
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dummy integration variable for z in equations (A37) and (A49) 


mR p 2 o 0 

defined in equation (A39) as 

m po n mdSe 

rR 2 a 0 

defined in equation (A33) as — 

m po I1 md®e 

distance between the fictitious compacted region determined by e d 
on the landing surface and that on the payload (or the payload 
itself, in the absence of penetration) ; evaluated in equation 
(5) 


variable used as d£ to define a length element of a rod of 
crushable material; used in developing equation (All) 


original mass of uncrushed crushable material; see figures 2(a), 
3(a), and 4(a) 


mass of material compacted by the landing surface; mass of V c i ; 
see figures 2 through 4 


mass of payload; see figures 2 through 4 

mass of material compacted by the payload; see figures 2(c), 3(c), 
4(c) 


first time derivative of mpi 

dimensionless mass defined in equation (A34) as 
dimensionless mass defined in equation (A44) as 


"Vo + m co - m cl 
m po 

m co _ m cl ” m pl 

m po 
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ma 
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n md 

n P 

n Pmax 
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q,q p 

q>q p 

q p 

\iax 5 ^Pmax 


dimensionless mass defined in equation (A42) as 


>0 + “pi 

m P o 


defined in equation (23) as 


defined in equation (24) as 


R P / Uo 2 \ 

R V * 2 e <J n Pn,axSe R P / 
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defined in equation (25) as — — — (see also eq. (26) 

cm g e^ 

for payload penetration with the simplified design) 
desired value of n p max 

maximum permissible value of n p max (2000 herein) 

g loading with iipg e as the payload deceleration 

maximum value of np that occurs during impact stroke 

absolute displacement of uncrushed crushable material after 
initial contact with landing surface; see figures 2 through 4 

first time derivatives of q and qp 

second time derivatives of q and qp 

absolute displacement of payload after initial contact between 
crushable material and landing surface; qp = q in absence 
of payload penetration; see figures 2 through 4 

values of q and q at end of impact stroke 


q s value of q (q = q p ) at which payload penetration occurs or 

would occur if unbonded 

R overall radius for spherical system 

Rp payload radius for spherical system 

r polar radial coordinate of point in A cl , where coordinate is 

measured in a horizontal cross section for spherical 
geometry; see figure 4(b) and equation (A25) 

rp polar radial coordinate of point in Ap X> where coordinate is 

measured in a horizontal cross section for spherical 
geometry; see figure 4(c) 
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SEA 


specific energy absorption; energy absorbed per unit weight of 
the absorber; defined in equation (8) in terms of the 
variables used herein 

dimensionless dummy integration variable used in equation (A36) ; 
= b 
" R 

dimensionless dummy integration variable used in equation (A45) ; 


t 
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V r 


v c ,v cl 

^max* Ci 
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cl 


V- 


PI 


v r 


time; used for time derivative interchangeably with dot 
velocity of the uncrushed crushable material; defined in 


equation (A16) as 


_ dq 
* = at 


max 


value of U at the instant the crushable material hits the 
landing surface (300 ft/sec herein) 

volume that would lie beneath the landing surface if there had 
been no crushing; shown in figures 3 and 4; defined in 
equation (A14) 

maximum values of V c and V c i , that is, the values reached at 
the end of the impact stroke 

compacted volume of material compacted by the landing surface; 


the volume of m, 


ci 


volume swept out by the payload as a result of the relative 
motion of penetration; defined in equation (A14) 

compacted volume of material compacted by the payload; the 
volume of mpi 

dimensionless velocity of the uncrushed crushable material; 

dq/dt _ U dy 


defined in equation (A54) as 


^ n md&e^p ^ n mdSe^- 


dx 


= v 


P 



value of v at the instant the crushable material hits the 
landing surface 


dimensionless velocity of the payload; defined in equa- 
tion fA541 as — — = 

dX dX dX 

value of v at the start of payload penetration, if any 
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original earth weight of uncrushed crushable material, m co g e 
earth weight of payload, m__,g„ 

pO 6 

dimensionless velocity of the uncrushed crushable material; 

defined in equation (A33) as — 

/n mdSe R 

value of w at the instant the crushable material hits the 
landing surface 

value of w at the start of payload penetration, if any 
dimensionless time; defined in equation (A39) as t 

dimensionless displacement defined in equation (A39) as 



dimensionless payload displacement defined in equation (A54) a s 
q p _ 

values of y and y p at end of impact stroke 
value of y at the start of payload penetration, if any 
dimensionless displacement defined in equation (A33) 
value of z at np = n p max 

value of z at end of impact stroke 
q P 

p max 

R 

value of z at the start of payload penetration, if any 

angle between normal to stressed area and direction of maximum 
normal crushing stress 

n Pmax _ ^ 
n des 


n md&e 

XT" 


e compacting strain of crushable material; more detailed 

definition given following equation (1) 
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fictitious value of assumed for design purposes; always less 

than e m 

compacting strain when it is uniform throughout hypothetical 

crushable material; determines the surface for stress evaluation 

and can be specialized in governing equations for simplifying 
approximations 

e m value of for actual crushable material (0.8 herein) 

0p angle shown in figure 4(c) and used in derivation of equation (A28) 

£ dummy integration variable for y 

p c k uniform density of hypothetical crushable material in swept-out 

volumes such as V c in figure 4; p c ]< can be specialized as zero 
to neglect variable mass without implying a massless crushable 
casing 

Pern uniform density of actual crushable material 

Pp payload packaging density; defined in equation (A43) as 

^po 

(4/ 3) irR p 3 

m po 

PpP defined in equation (A35) as 

(4/3) ttR 3 


a "mostly normal" and "mostly static" crushing stress of crushable 

material, averaged over maximum possible stroke prior to compacting 

a Q maximum of a as defined in equations (2) and (A24) 

Op 0 normal stress on the payload prior to penetration 

a pok a po under assumption Op 0 = Op^ = constant 

a v vertical component of o 

a v r dynamic value of a v as defined following equation (A7) 

a V p 0 vertical component of normal stress on the payload 

(p angle shown in figure 4(b) and used in derivation of equation (A27) 

4> angle shown in figure 4(c) and used in derivation of equation (A28) 
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APPROXIMATE ANALYTICAL MODELS FOR LANDING ENERGY ABSORPTION, 
INCLUDING THE EFFECT OF PENETRATION BY THE PAYLOAD 
INTO ITS CRUSHABLE CASING 
Robert W. Warner 
Ames Research Center 

SUMMARY 


Two approximate analytical models are defined for a landing configuration 
in which a spherical payload can sometimes penetrate into its crushable casing. 
Results for both models are found to agree reasonably well with two previous 
experimental measurements. Design examples are presented for an impact veloc- 
ity of 300 ft/sec. These are based on choices of zero or ’'perfect” payload 
bonding, and of either a balsa-like or honeycomb -like class of crushable mate- 
rial. The greatest difference between the models for these examples is a 29- 
percent discrepancy in the required maximum crushing stress. A particular 
pair of examples gives the unexpected result that penetration can provide a 
decrease in crushable material weight by a factor greater than 4 when the 
honeycomb-like class of material is required without penetration, but the more 
efficient balsa-like class is feasible with penetration. 


INTRODUCTION 


One means proposed for providing information on lunar and planetary 
surfaces consists of an unmanned instrumentation system that is hard landed 
(with or without terminal guidance) but designed to survive the impact. For 
such a landing, a crushable casing is one means for absorbing landing energy 
so the payload can survive and transmit information during and after the 
impact. The advantage of this approach depends greatly on how light the 
crushable casing can be made for a given impact velocity. 

Various aspects of the design of crushable casings have been 
investigated in references 1 through 9, but without including the effect of 
penetration by the payload into the casing. The primary purpose of the pres- 
ent paper is to evaluate that effect analytically. For this purpose, fairly 
general equations of motion are derived for the payload and the crushable 
material. These equations are then specialized for two approximate analytical 
models in which the payload and casing are spherical. The analytical models 
will be used in a number of design examples so that the effect of payload 
penetration can be evaluated and also one model compared with the other. The 
two models will also be compared with earlier models that do not permit 
penetration and with the test results of reference 10 and a private 

communication. 1 ___ 

^-Donald R. Cundall, December 1967. 
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OUTLINE OF THEORY 


Properties of Typical Crushable Material 

Material properties are an essential input to the theoretical 
development, A typical crushable material for landing impact energy absorp- 
tion has a stress-strain curve similar to curve ABCDE in figure 1(a). The 
material is elastic from point A to point B, and there is rebound between 
points D and E. (The rebound should be as small as possible.) Between points 
B and D there is a large volume change; and the crushing stress is relatively 
constant , with an average value a indicated by the horizontal dotted line 
in figure 1(a). Point C, where the crushing stress begins to rise abruptly, 
is called the compacting strain e and ranges from 0.6 for close-packed 
crushable materials to nearly 1.0 for open crushable structures. 

The area enclosed by the stress-strain curve ABCDE is the energy 
absorbed per unit volume of crushable material. This energy is maximized for 
a maximum permissible crushing stress (which determines the maximum landing 
vehicle deceleration for a given configuration) if the stress-strain curve 
approaches a rectangle. Hence the material is not strained beyond point D in 
figure 1(a), where the stress is equal to the prior maximum (point B) , even 
though the stress could go higher as shown. If the initial peak at B is too 
high to approximate a rectangle, it can sometimes be reduced by precrushing. 

A less desirable alternative is to accept the stress-strain curve but modify 
the load-deflection curve by changing the shape of the crushable material. 

Stress-strain and load-deflection curves for balsa wood, plastic foam, 
and honeycomb are given in references 11 through 14. These curves have been 
determined by dynamic crushing tests and by nearly static tests of specimens 
having uniform cross sections. Figure 1(b) is a sketch of a crushing test in 
which the specimen is compressed uniformly across the cross section by a 
plate, and figure 1(c) shows a test in which the specimen is penetrated by a 
plunger. In both tests, the material crushes in layers at a loaded surface, 
which may be at either end of the specimen in the case of the plate loading. 

It should be noted that a certain (exaggerated) amount of material is 
shown trailing outboard of the plunger in figure 1(c) due to shear resistance. 
This shear effect, as well as friction, causes a difference in the crushing 
loads of figures 1(b) and 1(c). In the case of balsa wood, however, refer- 
ence 7 indicates that this difference is small, while reference 13 indicates 
that it ranges from 5 to 15 percent (for a range of plunger sizes). 

It should also be noted that dynamic and nearly static crushing tests 
give different results. In fact, references 11 and 13 suggest ratios of 
static to dynamic crushing stress from 0.69 to 0.73 for various materials. 

Since the maximum velocity in the dynamic tests was 108 ft/sec, the dynamic 
effect is probably not due to variable mass (i.e., the accumulation of crushed 
material at a loaded surface) but rather due to damping and dynamic buckling 
(i.e., coupling between vertical and horizontal velocity) . Unfortunately, 
the effects of higher velocities on damping and dynamic buckling are not 
established for the present materials. 


2 



Summary of Fundamental Assumptions and Limitations 

In contrast to the typical material just discussed, the material 
assumed for theoretical analysis has a perfectly rectangular stress-strain 
curve (a so-called M rigid plastic” shape) . This curve is represented by the 
dashed lines in figure 1(a). It is shown bounded by the compacting strain e 
and the average crushing stress a of the typical material. In this case, 
the energy absorbed per unit volume will not match that of the typical material 
perfectly, but the boundaries of the rectangle can be adjusted slightly if 
need be. 

There are a number of other analytical assumptions that pertain to the 
crushable material, and there are several that do not. For convenience, all 
fundamental assumptions and limitations are listed in the present section as 
follows : 

1. Rebound is assumed absent. 

2. It is assumed that the effects of shear resistance, end fixity, and 
Poisson* s ratio are adequately incorporated in the analysis because 
of their presence in the crushing tests that determine the so-called 
"mostly normal’* and ’’mostly static” crushing stress (a, a Q , and a v ; 
see Notation) . 

3. The ’’mostly normal” and ’’mostly static” crushing stress a is also 
assumed to incorporate the dynamic effects of damping and dynamic 
buckling. This means that dynamic tests should be used to determine 
a (or static dynamic ratios, such as those deduced earlier from ref- 
erences 11 and 13, should be applied) . The velocities in the 
dynamic tests should not be large enough, however, to cause signifi- 
cant variable-mass effects, which are incorporated (when desired) in 
the present equations of motion. 

4. Shear deformations, such as the trailing of material outboard of the 
plunger in figure 1(c), are neglected. 

5. Separate vertical rods of material are assumed to crush vertically 
in the energy- absorbing process. 

6. The material is assumed to compress to the same compacted strain 
along any axis, regardless of the axis of maximum normal crushing 
stress . 

7. It is assumed that each particle in the uncrushed crushable material 
is moving at the same vertical velocity at a given instant and that 
the crushed material also has a uniform particle velocity, but a dif- 
ferent one from that of the uncrushed material. This implies the 
following corollary assumptions: 

a. Each successive layer of crushable material undergoes a 
jump in velocity as it moves from the uncrushed to the 
crushed region. 
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b. The elastic stress waves that established the uniform 
velocities must travel far faster than those velocities 
(i.e., the uniform velocities are subsonic). 

c. The deformation waves resulting from the elastic stress 
waves must be small enough not to affect the uniform 
velocities . 

8. For landing geometries in which there is doubt as to where the crush- 
ing by layers will start, it is assumed to start at the impacted end 
of the crushable material (in keeping with experimental observations, 
except at certain impact velocities too low to be of interest) . 

Thus, if figures 1(b) and 1(c) were considered to represent impact 
tests, the location shown for the crushed material requires that the 
material has been impacted by the plate or the plunger. If the 
crushable material had impacted the landing surface in figure 1(c), 
with the plunger resting on top of the material, there would have to 
be some crushing at the landing surface in conjunction with the 
plunger penetration (or M payload penetration 11 ) . 

9. It is assumed that there is no section of crushable material weak 
enough to permit penetration by any crushable material (i.e., pene- 
tration by anything but the payload) . The validity of this assump- 
tion is investigated in a later section. 

10. The density of the crushable material is assumed uniform in the 
uncrushed condition. 

11. The compacting strain of the crushable material is assumed uniform. 

12. Pure vertical translation is assumed. 

13. The analysis neglects all ringing and focusing of stress waves. 

14. A weightless exterior cover for the crushable material is assumed 
that is strong enough to prevent shattering of the crushable mate- 
rial. (For comparison with the experimental results of Cundall and 
of reference 10, however, the mass of the cover employed is assumed 
uniformly dispersed throughout the crushable material.) 

15. When there is payload penetration, that is, relative motion between 
the payload and the crushable material, the two are assumed perfectly 
unbonded. This rules out the interesting design possibility of pene- 
tration despite bonding and means there is no need to include the 
effect of tensile stresses over the upper surface of the payload. 

16. It is assumed that the landing surface is perfectly flat and 
perfectly rigid. 

17. Perfect rigidity is assumed for the payload. 
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Summary of Analytical Development 


The general vertically symmetrical landing geometry without tipover or 
horizontal velocity is illustrated in figures 2 and 3, which are used for 
developing the governing equations and for the definition of terms. (Note that 
with gravity terms being small, "vertical" can be any direction that is both 
parallel to the resultant impact velocity and normal to the landing surface.) 

In figure 2 the major limiting assumption is the absence of shear deformation 
(i.e., there is no compacted material dragged outboard of the payload and no 
compacted material lifted off the landing surface) . 

Evaluation of the stress integrals and variable masses in the governing 
equations is greatly facilitated by introducing the compressive compacting 
strain e and assuming that it is uniform throughout the crushable material, 
that is. 


e = = constant (1) 

Figure 3 shows the general vertically symmetrical landing geometry for zero 
shear deformation and the assumption of equation (1) . The latter assumption 
is illustrated by the fact that hpi is constant and m cl is a foreshort- 
ened image of V c . 

Figure 4 specializes the geometry to concentric spheres (although most 
of the resulting simplifications would be realized as well by concentric 
spherical segments). 1 2 It is further assumed that the crushable material has a 
specific direction for maximum normal crushing stress, that the material has 
been segmented and oriented to make that direction radial, and that stressed 
areas with normals differing by an angle a from that radial direction feel 
normal stresses a determined by the law 

a = a Q cos a , a < 90° (2) 

where the restriction a < 90° is required since the load has to be trans- 
ferred from the lower to the upper hemisphere of crushable material, and 
where the restriction is automatically met in all calculations for nonzero 
payload radius. (See ref. 5 for alternate anisotropic relationships.) 

Governing equations corresponding to the summary just given are 
developed in appendix A. They are specialized for various combinations of 
the following independent assumptions: 

1. Neglect of variable mass effect. In this assumption the accumulation 
of crushed casing material on the payload and/or the landing surface 
is neglected. The assumption is accomplished by setting = 0. 

2 The spheres were selected because of their ability to absorb impact 
energy in any direction (i.e., to handle unoriented impacts). Such impacts 
may occur due to aero dynamic ally unstable landing configurations, a strong 
laterial wind with a steep landing surface, or terminal guidance failure or 
absence. 
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Neglect of built-up material effect. In this assumption stresses are 
evaluated at the surface of infinitely thin sheets of crushed mate- 
rial rather than built-up volumes. The assumption is accomplished 
by setting = 1 . 

3. Neglect of variable resistance to payload penetration. In this 

assumption the stress and force on the payload retain their initial 
penetration value for the entire penetration stroke. The assumption 
is accomplished by setting Fp 0 (e) = 1. 


Appendix A is presented because a relatively complete theoretical 
development, including variable mass, built-up material, and payload penetra- 
tion, is needed in the literature. The development is relegated to an appen- 
dix because the details are not needed to understand the rest of the report. 
If the reader wishes to locate a specific result or derivation, he can refer 
to the various subdivisions of appendix A listed in the Table of Contents. 


DESIGN PROCEDURES FOR SPHERICAL GEOMETRY 


General Design Conditions 

The zero-velocity termination conditions defining the end of impact are 
equations (A46) and (A47) . For a minimum weight design, the termination con- 
ditions should occur when the payload or the compacted material built up on 
it (in the case of payload penetration) touches the compacted material built 
up on the landing surface, that is, when the sphere of radius Rp touches 
the region m C i in figure 4(b) or when the two compacted regions touch in 
figure 4(c). The size of the compacted regions is based on the true material 
compacting strain e m . The quantity c m is identical to the ck of equa- 
tion (1) except for the optional use of fictitious values for ek in the 
equations of motion of appendix A. If a margin of safety is desired, larger 
compacted regions can be envisioned on the basis of a fictitious design 
compacting strain called ed, where 


e d < £ m 


(3) 


If L is defined as the distance between the fictitious compacted 
regions at the end of the impact stroke, or between one such region and the 
payload, then the design condition for contact is 



where Rp is the payload radius and R the overall radius. With e^ 
replaced by it can be deduced from equations (A12) and figure 4 that 


(4) 


L = R - R p 


^Pmax 

e d 


(5) 
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where Rp max is the maximum absolute payload displacement during impact. 

if Qmax is the corresponding displacement of the uncrushed crushable mate- 
rial, then q p = q ma x i n the absence of payload penetration. When 

L > 0, the process is physically realizable, although not a minimum weight 
design for the present spherical geometry. When L < 0, the process is not 
physically realizable; but cases for L < 0 may be recorded in the process of 
seeking L = 0. 

A second design condition, having a less obvious relation to minimum 
weight design, is concerned with the maximum design deceleration n mc jg e , 
where g e is the acceleration due to gravity on earth and n m( j the maximum 
design g loading. If n p max £e t * le actua l payload maximum deceleration, 
the design condition is 


n Pmax < 
n md 


( 6 ) 


For certain types of energy absorbing material or structure, the minimum 
weight design calls for 

n Pmax = 
n md 


(7) 


For other types, any value of n p max / n md satisfying equation (6) may yield a 

minimum weight design, with the critical parameter being a property of the 
crushable casing (such as a Q ) . 

The quantity n P max / n md equations (6) and (7) is the maximum of the 

quantity np/n^j in equations (A32) and (A41) . (See also eqs . (A48) , (A51) , 

(A55) , and (A59) for various specializations.) If variable mass and built-up 
material are neglected in the analysis (by setting p c k = 0 and = 1, 
respectively) , then np occurs at q = q D = R/2 in the absence of payload 

penetration, providing q or qp becomes that large. 

A minimum weight design is sought for a given payload mass, mp Q , 
payload radius, Rp, impact velocity, U Q , design compacting strain, and max- 
imum permissible g loading, n m d, with the material choice being arbitrarily 
limited to two classes of crushable material. Where equation (7) gives the 
minimum weight design, the quantities to be determined are the overall radius, 
R, the maximum normal stress, a G , and the density, p cm , of the crushable mate- 
rial (from which the weight can be calculated) . If the optimization permits 
equation (6), then a property such as a G must be specified; and R and p cm 
remain to be determined (with n p max a by-product). In either case, the 

determination requires that a relation between a 0 and p cm be known for the 
crushable material, thereby effectively reducing the number of unknowns by 
one. 
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The required relation between a Q and p cm can be determined 
experimentally for a variety of materials and stated directly as in refer- 
ence 2 (p. 259). It is common, however, to introduce a parameter for which 
values are widely known, namely, the specific energy absorption (SEA) which 
is defined as 


/ °o \ 

SEA = e m ( (8 

m ypcmg e y 

and which should generally be as large as possible. The relation between 
o Q and p cm can be expressed by giving SEA in terms of o Q . With a = o Q 
and e = e m , the product a 0 e m is the area enclosed by the dashed lines in 
figure 1(a), which is the energy absorbed per unit volume. Dividing this 
product by p CIJ1 g e gives equation (8) and shows SEA to be the energy 
absorbed per unit mass. 

Figure 5 shows the variation of SEA with a Q for a balsa-like class 
of material and a honeycomb- like class of material (neglecting glue- joint 
weight and the effects of high impact velocity and low temperature) . The 
equation for SEA in terms of a Q given for the balsa-like material in fig- 
ure 5 is 

SEA = 24,000 ft- lb/ lb , 

SEA = 7,49 * 1C|6 ft- lb/lb , 

' 81 

The balsa- like material is balsa in the sense that the second of equations (9) 
is deduced from the curve fit of reference 2 (p. 259) which closely approxi- 
mates the data therein for soiled balsa of various densities. The material 
is considered to be only balsa-like, however, because the first of equations 
(9) and the dashed portion of the balsa-like curve in figure 5 are assumed to 
be valid for a hypothetical cored balsa (for which cores of material are 
removed in the radial direction) . The assumption is that roughly one-third 
of the material can be removed from the lightest (lowest density) solid balsa, 
giving a o Q range from 800 to 1,200 psi, without reducing the SEA by 
introducing significant buckling, end effects, and/or Poisson's ratio effects; 
and this assumption seems reasonable. It should be noted that the only cor- 
roborative case where SEA decreases with increasing a Q in reference 4 is 
for dry balsa (0-percent moisture), atmospheric pressure, and an ambient tem- 
perature of approximately 78° F. 

The curve for the honeycomb- like material is presented in figure 5 so 
that results can be deduced for lower SEA values and also for the relatively 
common case where SEA increases with increasing a Q . The equation given in 
figure 5 is 

SEA = 478.5 a 0 0 ' 446 ft-lb/lb , 600 psi ^ a Q $ 1700 psi (10) 


800 psi = a Q = 1,200 psi 


( 9 ) 


1,200 psi = a 0 1 1,800 psi 
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The material is called honey comb -like rather than a specified honeycomb 
because equation (10) is deduced from a curve fit in reference 2 (p. 259) 
that rather loosely fits data for several types of aluminum and fiberglass 
honeycomb having a variety of densities. It should be noted that reference 2 
makes no mention as to whether equations (9) and (10) incorporate the dynamic 
effects of damping and dynamic buckling (see item 3 in "Summary of 
Fundamental Assumptions and Limitations") . 


Design Procedure for Simplified Model Without Payload Penetration 


For the analytical model labeled "simplified" in the present report, it 
is assumed that variable mass can be neglected (setting p c k = 0) , that built- 
up material can also be neglected (e^ = 1) , that gravity forces are insignifi- 
cant (g = 0) , and that the resistance to payload penetration is constant 
( F po(e) = 1) • Under the first three of these assumptions (with the last 
required only for penetration) , equations (A55) and (A62) apply (in a simpli- 
fied form with g = 0) and the division of the latter by the former gives 
(with definitions from eq. (A54)) 


v r 


U 2 
u o 


max 


n p/ n md n p& e ^p 


CR/R p ) - [C2/3)y max ] 

y [ (R/Rp) - y] 


(ii) 


Equation (A55) is a parabola with n p max at y = (R/Rp)/2. Hence, for the 
present restrictive case, n P max Se t ^ ie acceleration at as long as 


max 


y max ~ (R/Rp)/2. If y max ^ (R/R p )/2, n pmax g e is the acceleration at 
y = (R/Rp)/2. Thus equation (11) can be expressed in terms of n p max as 


U 0 2 




W [1 " (2/3)(R p /R)y max ] 


1 “ ( R p/ R )ymax 


max 




Rr 


U 2 

5 4 — y 2 l 1 

n Pmax g e R p R max V 


- — y ) 

3 R y max/ ’ 


max 


1 L 

2 R r 


(13) 


If equations (12) and (13) are multiplied through by Rp/R and if the first 
'' ' (R/Rp)) the result is 


of equations (A54) is used (giving y = z 

max max 




z [1 - (2/3) z ] 
max L v ' max J 


n Pmax^e R 


1 - z 


* I 


max 


max 


(14) 
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- X , -«* 2 I 

n Pmax s e R " ax| 



z 


max 



max 


1 

2 


(15) 


Equations (14) and (15) , which are applicable only in the absence of 
payload penetration, are plotted in figure 6 as z max versus U 0 2 /(np max g e )R. 


Figure 6 can be used to determine 
and n p max 5 and the payload radius 


z max = %ax/ R on the basis of u o> R > 

Rp simply has to be small enough not to 


interfere with 


q 

hnax 


An interesting feature in figure 6 is that z max I s determined with no 
knowledge of cr 0 , m co , or mp 0 . Once z max is known, however, m co can be 
found for a given mp Q and o Q according to 


m 


co 


2 7 TO R 3 

Q 2 I 

2 Z max 

lb/ 


1 - — z 
1 3 nnax 


- m 


po 


(16) 


as derived from equations (A62) and (A39) with g = 0. It is still not 
necessary to know Rp. For a minimum weight design, however, Rp must be 
known or determined; and equation (4) must be satisfied according to the 
definition in equation (5). In addition, the SEA , as defined in equation 
(8), may be known instead of a Q . When equations (4), (5), and (8), with 
z max 5 q max/ R = q Pmax^ R and mc0 = 0/3)irp cm (R 3 - R p 3 ) , are introduced into 
equation (16) , the result is 


J 


ma 


g e (SEA) 




1 - (z / e ,) + [(z / e ,) 2 /3l 

max rr L mar rr / 


max 


z [(1/2) 
max L ^ J 


(z / 3) ] 

^ max 


(17) 


where the symbol J mo is introduced for convenience. Equation (17) can be 
used to determine m co for a given z max (which implies knowledge of U G and 

g e as indicated in fig. 6) if ed, e m , SEA, and mp Q are known (with e m 

actually canceling the same quantity in SEA) . In figure 7, equation (17) is 

plotted for ej = 0.7, 0.8, and 0.9. The three cross-plotted values of Rp/R 

serve as a reminder that figure 7 represents a design for which contact would 
occur between the payload and the compacted material if that material had a 
compacting strain of ej instead of e m (see the sketch, fig. 7) . 

In the absence of payload penetration, figures 6 and 7 are sufficient 
for the simplified model if R is given and Rp is to be determined. If 
Rp is given, figures 6 and 7 remain useful as a check and as a means of 
determining z m ax* but two different figures, based on Rp, are more useful 
for the original design. The first of these is determined by introducing 
equations (4) and (5) into equations (13) and (12) as a substitute for 
y m ax E q max/ R p = q p max / R p> yielding 
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U o 2 4e d 2 [l - (R p /R)] 2 jl - (2/3)e d [l - (R p /R)]j 


n_ g R 
Pmax e P 


R p /R 



1 

2E d 

(18) 


U 


n n g R 

p max e p 


M 1 


(Rp/R)] ( 1 - (2/3) ^[1 - (R p /R)] 
(Rp/R) j 1 - s d [l - (Rp/R)] 




(19) 


Equations (18) and (19) are plotted for ej = 0.7, 0.8, and 0.9 in figure 8. 

Figure 9 is the companion to figure 8 and is found by using equations 
(4) and (5) to define a replacement for z max E ^max/^ = ^Pmax^ in e 0 ua ~ 

tion (17). The Rp/R found in figure 8 then determines J ma in figure 9, 
and J ma determines m co as before. The sketches in figures 8 and 9 show 
the same design configuration as that in figure 7, leaving figure 6 as the 
only one for which contact is not required for compacted material based 
on 


Design Procedure for Simplified Model With Payload Penetration 

Payload penetration was not permitted in figures 6 through 9 (because 
of a hypothetical perfect bond). If penetration is now permitted (because of 
the total absence of bonding) , the first step is to determine whether it will 
occur. It will occur if the design without penetration gives 


z > 
max 


( 20 ) 


where z s is the dimensionless displacement at which the area of material 
being crushed becomes large enough to cause sufficient deceleration for pene- 
tration. The quantity z s is defined by equation (A58) as modified by the 
change of variable, y s = z s (R/Rp) . The modified equation, plotted in figure 
10, is 



+ 


m co\ 

m po J 


2z s (l 


- z s) 


( 21 ) 


The curve is cut off at z s = 0.5, which is the maximum value at which pene- 
tration can occur for the present approximation (as seen by the maximum force 
in eq. (A52)). Figure 10, while not ultimately essential for design purposes, 
is useful as a check, and the magnitude of z s is significant in evaluating 
the importance of penetration or potential penetration. 
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The first of the two main design figures for penetration is based partly 
on the first of equations (A59) and the last of equations (A39) with g = 0 
and n p = n Pmax = constant. This yields, as expected for the present 

assumptions , 


n 


*y 


Pmax mpoge 


(22) 


Equation (22) is used, together with equation (21), equation (A64) for g = 0, 
the last of equations (A39) , the first and second of equations (A54) , and 
equations (4) and (5) for z Pmax = Ppma)/ 1 *’ to identify the variable Nru as 
follows: 


N 



L u ° 2 )-i z * 

1 - (4/3) z s 

l 2e .n^ g R e, 

\ d p max e y / d 

_2(1 - z s ) _ 


(23) 


The purpose in isolating Nru is to obtain a single unknown (R p /R) in terms 
of z s and ej, with a factor that is a function of known quantities (U 0 , e^, 
n Pmax’ ® e ’ **p) • T ^ ie next ste P is to isolate a variable for m co /mp 0 without 

R p /R. This is accomplished by dividing equation (21) by the square of equa- 
tion (23) . The resulting variable, called N mu , is 

1 + (m c 
N 


‘co/ m po^ 


2z s d - 


z s^ 


mu 


[l + W 0 2 /2c d n Pmax g e R p )] 2 ^ _ (Zs / Ed) J tl _ ( 4 /3 )z s ]/ 2(1 - z s )})' 


Figure 11 is a plot of Nru versus N mu for = 0.7, 0.8, and 0.9. 

The plot is constructed by selecting numbers for z s between 0 and 0.5 and 
calculating the corresponding values of Nru and N mu according to equations 
(23) and (24). A relationship between R and m c0 is established in figure 
11 in terms of the known quantities U 0 , cd> n p > g e > R p> and m po (i- n con ~ 

trast to fig. 8 without penetration, where R can be determined from known 
quantities and used in determining m co in fig. 9) . The sketch in figure 11 
shows the design condition implicit in equations (4) and (5) for penetration, 
namely, contact between the two regions of compacted material if e m were 
replaced by ej. 


A definition and a volume density relationship are useful at this point, 
namely: 


N 

ma 


Wpo ( m co/ m po) 

tfcmM’V 


w 


co 


m 


co 


(p cm g e ),,,< p 3 W 3 



1 


(25) 


where W p0 and W c0 are the original weights of the payload and crushable 
material, respectively. Equation (25) is plotted in figure 12; it is useful 
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with or without payload penetration and does not depend on any assumptions . 
For the case of payload penetration, however, and for the special assumptions 
of the simplified design, equation (22) applies and can be substituted into 
equation (25) with equation (8) to give 


= SEA / mco \ 
n Pmax £ m R P \ m P ° J 


(26) 


With equation (26), m C0 /mp 0 can be determined from N ma without knowing Pern- 
Thus figure 12 becomes the companion for figure 11 in an iterative design pro- 
cedure for payload penetration. The procedure is simply to select an initial 
value for m co /mp 0 , calculate N mu in terms of ni C0 /mp 0 and known quantities, 
determine Nru from figure 11, evaluate Rp/R in terms of N#jj and known 
quantities, determine N ma from figure 12, calculate a second value of 
m C o/ m po in terms of N ma and known quantities, and repeat the process until 
two values of m co /nip 0 agree to the accuracy permitted by the figures . 

The iterative procedure just described could have been avoided, of 
course, by combining equations (21) through (26) into a polynomial for z s . 

This polynomial could be solved for selected parametric values of 
U Q 2 / 2 e d n p max g e Rp and SEA / n P max e m R p ^ together with values of cj such as those 

selected for the iterative charts, and Rp/R and m C0 /mp 0 could be determined 
accordingly. Such a procedure was avoided, however, because of the strong 
likelihood that a reasonably limited group of parametric values would not have 
sufficiently broad applicability. 


Specialization of Design Procedures for Simplified Model 
to Materials of Figure 5 

Regardless of whether penetration is absent or present, it is apparent 
that all of the design figures (6 through 12) are essentially independent of 
the material or structure selected for energy absorption. To calculate m c0 
from figures 7, 9, and 12, however, an SEA has to be selected. With m co 
and R determined by the figures and Rp known in advance, p cm can be calcu- 
lated (with the aid of fig. 12 if desired); and o Q can then be found for the 
selected SEA according to equation (8) . All this implies the assumption 
that a material having the calculated properties is available or that a 
corresponding structure can be constructed. 


For the comparison purposes of this paper, however, the materials are 
restricted to a choice between the two classes described in figure 5. With 
figures 6 through 12 having been constructed independently of figure 5 (in 
the interest of generality) , the use of figure 5 imposes a trial-and-error (or 
transcendental) solution for those designs in which payload penetration is 
absent. The trail-and-error solution is aided by incorporating equations (8) 
and (25) into equation (17) to yield 


SEA 


e m J ma 

(ge^d/Uo 2 ) - (Jma W po/°o N ma^p 3 ) 


(27) 
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On the basis of given values of e m , g e > e d, u o> w po> R p> and n p max > P lus 

values of J ma and N ma determined according to figures 8, 9, and 12, the 
quantities e m J m a , g e e d/ u o 2 > and J maW po /Nina , n' R p 3 in equation (27) can be 
calculated in advance. Equation (27) and figure 5 then become the basis for 
the numerical solution. 

A sample calculation for the trial-and-error solution just described is 
given in appendix B. It should be noted that the trail and error of appendix 
B would have been eliminated if the SEA had been constant - despite varia- 
tions in a 0 , as for the balsa-like curve in figure 5 when a Q < 1,2 00 psi. 
There is, however, another type of trial and error that occurs when a spe- 
cific value of a G , rather than n p max > i s sought. In this case, calcula- 
tions like those in appendix B (but without the o Q - SEA trials) must be 
performed for successive selections of n p until o 0 approaches the 
desired value. maX 

If (in contrast to appendix B) payload penetration is present, the use of 
the specific materials in figure 5 does not impose a trail-and-error solu- 
tion. The reason is that the unknown m co does not appear in equation (22) . 
In addition, the simplicity of equation (22) means that either o Q or n p 
can be selected without requiring trial and error. 

A sample calculation for the iterative penetration procedure described 
earlier is given in appendix C for the balsa-like material of figure 5. 

Only three iterations are needed in appendix C because of a fortunate initial 
guess of m co /m po . The initial guess was based for all calculations on other 
penetration cases or prior examples without penetration, and the worst guess 
required five iterations . When a penetration design is not feasible, either 
the o 0 value of equation (22) will be beyond the range of figure 5, or the 
iterations will move off the curves of figures 11 and 12. 


Computer Procedures 

A set of three computer procedures has been programmed to evaluate the 
governing equations of the impact problem for spheres (eqs . (A32) , (A37) , and 
(A41)) in accordance with the appropriate termination conditions (eqs. (A46) 
and (A47)) and the appropriate design conditions (eqs. (4), (6), and (7)). 
These procedures are described in appendix D. They can be used to check the 
simplified model designs or to initiate simplified designs and other designs 
based on more complicated models. 

Specifically, of the three computer procedures, two are designs in that 
they automatically iterate initial guesses to determine required crushable 
casing parameters. One of the two design procedures varies the overall radius 
R and the material maximum crushing stress o Q to achieve a desired accel- 
eration, as indicated by equation (6) or equation (7), and a desired ratio of 
stroke to available stroke, as suggested by equation (4) . This program is 
applicable only in the absence of payload penetration. The other design 
varies only R for a selected a 0 (fixed material once SEA or a plot of 


14 



SEA versus o Q is selected) and achieves only the desired stroke ratio. It 
applies both with and without penetration. 

Penetration requires only a search for R since o Q is essentially 
determined in that case regardless of R or the stroke. The reason becomes 
apparent when the first of equations (A41) is modified at the start of pene- 
tration (e = de/dx = 0) by equations (A45) , (A42 ) 9 and the last of equations 
(A39) to give the expected result 

irRp 2 CT 0 

n p m + (2 

Pma x g e m po g e 


where the approximate equality sign is used because n p max differs slightly 

from the value at the start of penetration for the detailed model (by at most 
one part in a thousand for the present examples) . Equation (28) determines 


for a given n T 


max 


regardless of R or the stroke. 


The third computer procedure can also be used with or without penetration 
but is not programmed to iterate and produce a design. Hence it can be used 
only to check a given configuration. 


It should be noted that all three programs permit SEA to be selected 
arbitrarily or calculated (after a Q is selected) according to equations (9) 
and (10) for the materials considered herein. 


Design Procedure for Detailed Model With and 
Without Payload Penetration 


The detailed model is the second of the two approximate analytical 
models. Where the simplified model had p c ^ =0, = 1, g = 0, and 


Fp 0 (e) " 1 j the detailed model has 
equation (A45) , 

F fe) = 
po 


Pck = P 


cm* 


e k 


e m> g = §L> and from 


/ 


(eSp 

/2es T 


l)s n ds T 


+ e + 1 


The Pck = Pern equation means that the detailed model incorporates variable 
mass (accumulating on the payload and/or the landing surface) according to 
the crushable material density, and = e m implies a finite volume of 
compacted material instead of an infinitely thin sheet to determine the sur- 
face for stress evaluation. The equation g = gL simply incorporates a 
negligibly small gravity term for completeness, and the integral for Fp 0 (e) 
permits a calculated deviation from a constant resistance to penetration. 


While the computer procedures described in the previous subsection 
constitute an alternative to figures 5 through 12 for the simplified model, 
they constitute the only design procedure presented herein for the detailed 
model. Of course, figures 5 through 12 remain useful as starting points for 
the computer iterations . 
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DESIGN RESULTS AND DISCUSSION FOR SIMPLIFIED AND DETAILED 
MODELS HAVING SPHERICAL GEOMETRY 


Description of Design Examples and Most General Results 

Thirteen examples have been calculated by the simplified and detailed 
analytical models described in the previous section. The impact velocity 
U 0 is taken to be 300 ft/sec for all examples, the material compacting 
strain e m to be 0.8, the design compacting strain to be 0.7, and the 

payload maximum g loading n p max to be 2000 or less. The examples vary 

in the given payload weight Wp 0 , the given payload radius Rp, the presence 
or absence of a hypothetical perfect bonding, and/or the given material. 

All examples are based on a choice between two classes of crushable 
material - the balsa- like and honeycomb- like materials defined in equations 

(9) and (10) and shown in figure 5. The desired value of n D (for which 

*max 

the label n^gg is used where needed) is given for some examples and the 
maximum crushing stress o 0 for others. 

The results for all examples are presented in table 1 in terms of the 

presence or absence of penetration, the value of n n (if not given) or of 

p max 

a G (if not given), the crushable material density P C m§e> the specific energy 

absorption SEA, the stroke-to-potential-start-of-penetration ratio q D /q s , 

Fmax 

the overall radius R, the crushable material weight W co , the total weight 
W co + Wp 0 , and the dimensionless unused stroke L/R. When L/R is negative, 
the payload would go too far (slightly, for the present examples) and cause 
excessive accelerations if it were not for the stroke margin of safety given 
by the use (in defining L/R) of the fictitious design compacting strain 
rather than the material value e m . 

The resulting quantities just listed are presented for the simplified 
and detailed models in table 1; the corresponding ratios of detailed to sim- 
plified results are presented when the results are numerical. It is apparent 
from the W co ratios that the simplified model has the lower W co . 

An important ground rule for the ratios of table 1 is that either the 

a D ratio or the n D ratio is required to be essentially unity. The 

*max 

choice is made, after the simplified model has been calculated, in favor of 
the lightest resulting detailed model. This gives the best weight ratio, 
that is, the ratio closest to unity. 

In view of this choice, it is not surprising that the worst o 0 or n n 

p max 

ratios (of detailed to simplified results) are farther from unity than the 
worst weight ratios. In fact, the worst ratio of all is 1.2917 for o Q in 
case 1 (except for L/R, where numbers approach zero, fortunately) . This 
ratio, although not excessively different from unity, represents a large 
enough difference in required material (within the honeycomb- like category of 
case 1) that a designer would wish to know whether the simplified model is 
numerically more realistic than the detailed model or vice versa. 
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In a footnote of table 1 it is noted that the simplified results are 
"chart 11 results (derived from figs. 5 through 12) except for the L/R column 
and cases 7 and 12, which are automatic computer results. The two cases 
(both without penetration) were added after the computer design programs 
became available and were calculated by those programs for convenience. 

The simplified results determined by the chart method have been checked 
by the computer checking program (as opposed to a design program) . The 

derived computer quantities n D and q D /q s had a maximum error of 1 

F max F max 

percent relative to the chart method, and the derived computer quantity L/R 
was close to zero (ranging from -0.001781 to 0.001163 with an average of 
-0.000194) as compared to a zero chart value. Slightly larger differences 
between the chart method and the computer checking method occurred for W co 
and W CQ + Wp 0 , but these differences resulted purely from errors in calcu- 
lating crushable material volumes in the chart method. The chart and compu- 
ter results combined to form smooth plots for the simplified model, as seen 
in figures 13 and 14. 

All numerical results from table 1 are plotted in figures 13 and 14 
except for W co + Wp Q (which is considered less important than W co ) , L/R 
(which is often sporadically variable in sign), cases 3, 5, and the detailed 
model for case 2 (which will be discussed later) . Note that the curves for 
the simplified and detailed models are roughly parallel (and only slightly 
curved) over their mutual abscissa range. Hence trends are the same for the 

two models. Even the crossover in figure 14(d) for q^ /a without pay- 

Pmax * 

load penetration is almost parallel and thus maintains the trend with only a 
slight comparison reversal between models. 


Effect of Payload Radius (and Payload Packaging Density) 
for Balsa- like and Honeycomb- like Materials 
Without Payload Penetration 

Crushable casing properties and performance are presented in figure 13 
(with and without penetration) as functions of payload radius Rp (with four 
numbers attached for payload packaging density Ppg e ) for honeycomb- like 

material, a payload weight of 100 lb, and an n D value of 2000. The cor- 

* max 

responding cases in table 1 are 1, 4, 6, and 7 without penetration and case 2 
(simplified model only) with penetration. 

In the absence of payload penetration, figures 13(a), 13(b), and 13(c) 

show that an increasing Rp (decreasing Ppg e ) causes W co , SEA, R, a 0 , and 

PcmSe to increase (with n D held at the maximum permissible value of 

r max 

2000, which minimizes W co according to preliminary calculations). The poten- 
tial penetration ratio qp^/q s exists only for the simplified model at 

Rp = 0.6 ft (among the plotted points), indicating penetration to be impossible 
for the other cases. Between the two models, the simplified model has the 
lower W co (slightly), the lower SEA, the higher R, the lower o 0 , and the 
lower PcmBe- The lower p cm g e is the only apparent reason for the lower 
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W co . It results from the lower o Q , which is made possible by the fact that 
the simplified model sees a higher average stress for a given cr 0 than does 
the detailed model. 


Somewhat similar curves are shown in figure 14 for balsa-like material 
and a payload weight of 450 lb (the cases from table 1 being 8, 10, 12, and 13 
without penetration, and 9 and 11 with penetration). It should be noted that 
the 450- lb payload, in contrast to the lighter payload, permits the use of the 
efficient balsa-like material without causing excessive decelerations in the 
absence of penetration. When penetration is absent, results are presented for 
the lowest o 0 considered, namely, 800 psi. According to preliminary calcula- 
tions, this a Q value gives minimum weight results for the balsa- like material 

(cored when o n is less than 1200 psi) , as opposed to n D = 2000 for the 

F max 


honeycomb- like material. The stress of 800 psi determines SEA as 24,000 
ft- lb/lb and P C mSe as 3.84 lb/ft 3 . Figures 14(a), 14(b), 14(c), and 14(d) 
then show that an increasing Rp (decreasing Ppge) causes W co , R, and 

" ^ q_. /q_ to decrease (decreasing the likeli- 

Fmax b 


n D to increase but causes 
F max 

hood of penetration, as expected for increasing Rp) . Note that the 

q D /Qs of 0.730 at R_ = 1.6 ft is the only case where penetration is 
F max F 

impossible because q-n /q_ is less than one. A comparison of the simpli- 
fied and detailed models indicates that the former has the lower W co , the 


lower R, and the higher n-r, . The q- /q c curves have a shallow cross- 

F max Pmax 5 

over (in the absence of payload penetration) , and q /q,, does not exist 


for the detailed model at 


Rp = 1.6 ft. 


Pmax 


This absence of 


/q s for the detailed model when it exists for the 


q Pmax' n s 

simplified model has been a recurring theme without penetration, as seen in 
table 1, and indicates that the detailed model is the least susceptible to 
penetration. The reason for this is the relatively low maximum g loadings, 


n 


Pmax ’ 


for the detailed model in figure 14(c). 


Effect of Payload Penetration 

The data for payload penetration in figures 13 and 14 are limited by the 

fact that penetration often does not occur even when the payload is unbonded 

(see the q n /q_ column in table 1 for case 2 with the detailed model and 
F max 5 

for cases 4, 6, 7, and 13 with both models) . It is seen in figure 13 for 
honeycomb- like material, that penetration reduces W CQ , increases SEA and a 0 , 
decreases R, increases Pcm?e> and decreases qp /q s for the simplified 

model at Rp = 0.6 ft. The impossibility of penetration for the plotted 
abscissas higher than Rp = 0.6 ft means that the W C q curve for the higher 
Rp values could have been combined with the penetration point at Rp = 0.6 ft 
to form a W c0 curve for an unbonded payload, with an obviously beneficial 
effect. 

The W co benefit due to penetration at Rp = 0.6 ft results from the 
increase in SEA and a D and the corresponding decrease in R. These 


18 



quantities can change because the honeycomb- like material is allowed to change 

within its class to maintain n-n = 2000 for penetration (which gives the 

F max 

lowest W co , according to preliminary calculations, as it did without 
penetration) . 


In figures 14(a) and 14(d) for balsa-like material, it is apparent that 
payload penetration increases W co and decreases qp /q s for both models 

at Rp - 1.0 ft and R p = 1.2 ft when SEA, a G , and p cm g e are constant at 

24,000 ft- lb /lb , 1,200 psi, and 5.76 lb/ft 3 , respectively (where o Q = 1,200 
psi defines the elbow of the SEA curve in fig. 5 and gives the minimum W co 
according to preliminary calculations) . The overall radius R is shown in 


figure 14(b) to be increased by penetration at Rp = 1.0 ft and decreased at 
Rp = 1.2 ft; the apparent contradiction between tne R and W co effects of 

-L J. ■ * _ -I J T t 1 * ~ 1 - 


penetration is resolved by recalling that the lowest W co without penetration 

was for a lower PcmSe* namely 3.84 lb/ft 3 , than for penetration. In figure 

14(c) the effect of penetration on n n is seen to be the same for the two 

u max 

Rp = 1.2 ft. 


models at R r 


Pmax 

= 1.0 ft but different at 


The W co effect of penetration indicates that a perfect payload bonding 
would be desirable for the cored balsa- like material at Rp = 1.0 ft and 
Rp = 1.2 ft. In the event that such a bonding is not feasible or trustworthy, 
however, the case of an unbonded payload must be considered. Hence undefined, 
wavy- line transitions between penetration and no penetration are shown in 
figures 14(a), 14(b), and 14(c), but not in figure 14(d) (due to lack of 
space). The most important of these transitions is for W co in figure 14(a). 
If the W co transition were specified in the area of the wavy line, it would 
define a curve for an unbonded payload since the wavy line skips over 
Rp = 1.4 ft, where bonding is required to prevent penetration at 
o 0 = 800 psi. 


Stresses between 800 and 1200 psi (the minimum weight value for 
penetration) will presumably be useful in the transition. Thus a variety of 
specified transitions will be possible. For an unbonded payload, it would be 
desirable to seek a minimum W co in the transition region for a single o 0 
at which penetration is ready to begin at the end of the stroke; but this is 
beyond the scope of the present report. The reasonable assumption, however, 
is that such a minimum W co exists and indicates an important design trade- 
off for cored balsa, a point where a further increase in payload packaging 
density (the abscissa of fig. 14(a)) is undesirable. 


The most important effects of payload penetration discussed so far have 
been a decrease of W c0 for the honeycomb-like material and an increase of 
W co f° r the balsa-like material. Both comparisons have been based on mate- 
rial variations for minimum weight within the categories. If these varia- 
tions are not allowed, penetration can always be expected to increase W co . 

Except for the roughly 7-percent decrease in W co for the honeycomb - 
like material at Rp = 0.6 ft, there has been no advantage of penetration 
reported up to this point. If, however, a design is permitted to change 
from the honeycomb -like category to the balsa- like category, then a major 
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penetration advantage occurs in the present examples for the 100-lb payload. 
Specifically, a comparison of case 1 (Rp = 0.6 ft with n p max = 2000, 

honeycomb- like, no penetration) and case 3. (Rp =0.6 ft, a Q = 1200 psi, 

balsa-like, penetration) from table 1 shows a weight saving due to penetra- 
tion by a factor greater than 3 in W co and greater than 2 in W co + Wp Q . 

In addition, a comparison of case 4 (Rp = 0.7 ft with n p max = 2000, 

honeycomb- like, no penetration) and case 5 (Rp =0.7 ft, a 0 = 902.4 psi at 

n Pmax 2000, balsa-like, penetration) shows a weight saving due to 

penetration by a factor greater than 4 in W co and almost 3 in W c0 + Wp 0 . 
Presumably, a still more impressive comparison would have occurred at 
Rp = 0.8 ft except that a penetration design could not be achieved for 

n Pmax * 2000 • 

The honeycomb- like material is used in the comparisons when penetration 
is absent. The reason is that the balsa-like material, having the higher 
SEA over the same stress range (fig. 5), is considerably lighter than the 
honeycomb- like material. In fact, with the force determined by the stress 

range, the balsa-like material provides so little mass that the g loading 

exceeds the limit of 2000 for the 100- lb payload; and the heavier honeycomb- 
like material is required in the absence of payload penetration. (Note that 
this is not true for the 450- lb payload.) 

When penetration is permitted, however, for the balsa design by removal 
of the hypothetical perfect bonding between the top of the payload and the 
crushable material, then the only stresses acting on the payload are the 
crushing stresses at the bottom of the payload. This reduces the force suf- 
ficiently that the payload maximum g loading can be held to 2000 even for 
the balsa-like material with the 100-lb payload. The high SEA of the balsa- 
like material then produces the large weight saving due to penetration. 

Note that this weight saving applies for the selected classes of 
material; it would obviously be decreased if the honeycomb -like material were 

replaced by an intermediate class just heavy enough to bring n n down to 

F max 

2000 without penetration. The possible availability of such a class is 
indicated by the less efficient balsa reported in references 4 and 5. These 
balsas may be made still less efficient, in the sense of being heavier, by 
glue joints or by the addition of weights (a required weight being the cover 
for the crushable material) . 

Thus the large weight saving due to penetration is clearly restricted 
to the presently selected materials. It is obviously, however, a phenomenon 
worth considering. The fact that a heavy honeycomb material is still being 
considered in recent design studies is indicated in references 5, 8, and 9; 
and the g loading issue raised herein is emphasized in reference 5. 


Comparison With Previous Analytical Models 

The foregoing results apply to the specific analytical models 
considered here; and a question arises as to their validity. They can be 
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considered valid in one sense if they constitute a logical analytical 
extension of a reasonably standard body of theory having at least a limited 
experimental verification. The extended analytical models can then be compared 
with future experiments to establish experience factors as well as means of 
improving the models. 

The logical extension of standard theory is considered first, with the 
limited experimental verification to follow. As pointed out in the Introduc- 
tion, the present analytical models extend a significant amount of prior work 
(e.g., refs. 1 through 9) to include the effect of payload penetration. The 
summary of assumptions and limitations given under "Outline of Theory" is the 
same as the assumptions in references 1 through 9, where stated, except for a 
few variations on the assumption of equation (2). 3 

With the fundamental assumptions the same, it is no surprise that the 
basic equation of motion (with the anisotropy relation removed) is the same 
for the present simplified model without penetration as for comparable models 
(mass assumed constant, infinitely thin disk of crushed material, gravity 
neglected) in references 1 through 9, where stated. For the present detailed 
model, the most nearly comparable and completely described prior model is 
given in reference 2. The two models are identical except that the model of 
reference 2 has a uniform and isotropic crushing stress instead of the par- 
ticular radial distribution of equation (2). In appendix E, it is shown that 
such a crushing stress converts the fundamental equation (A19) into equation 
(E7) , which agrees exactly with the fundamental equation (1-4) in appendix A 
of reference 2. 

The analytical extension to include payload penetration is logical for 
the detailed model in the sense that no fundamental assumptions are added to 
the list given under "Outline of Theory." For the simplified model, the 
assumption is added that the force resisting penetration is constant, but it 
has already been pointed out that this is true for the detailed model to one 
part in a thousand or less, for the examples considered herein. 


Comparison With Previous Experiment 

The applicable experimental information known to the author is limited 
to two configurations, one tested at an impact velocity of 374 ft/sec (ref. 

10) and the other at roughly 220 ft/sec (footnote 1) . The latter velocity 
actually represents an average of velocities ranging from 215 to 225 ft/sec 
for four nearly identical tests. At the higher impact velocity, only a 
single test is considered because the exterior cover for the model in that 
test was the only one (out of four) sufficiently strong and resilient to 

prevent model disintegration. __ 

3These variations in anisotropy can be important for various crushable 
materials, as pointed out in reference 5. It may be desirable to consider the 
standard theory as incorporating the variations of reference 5, with the 
assumption of equation (2) considered as an example for evaluating the effect 
of penetration. If, however, it is desired to have only one anisotropy rela- 
tion for a standard theory, then the assumption of equation (2) is desirable 
since it is simple and widely used. 
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The tests just described demonstrated the existence of payload 
penetration, depending on the efficiency of the payload bonding. They also 
indicated that both the detailed and simplified models provide reasonable 
(within 10 to 30 percent) estimates of measured impact deceleration and 
stroke. There are not enough data to determine which model is quantitatively 
better. 


Alternate Models 

The analytical models considered here are manifestly only two out of 
many possibilities. Even with attention restricted to the present computer 
programs, it would be possible to investigate separately the effects of the 
limitations calling for p c ^ = 0 and efc = 1, rather than both together (the 
assumptions of g = 0 and Fp 0 (e) = 1 being trival) . It would also be possi- 
ble to set pj-k = 0 in the present programs for integration of the equation of 
motion but not for evaluation of the acceleration. The programs could be mod- 
ified to incorporate various levels of bonding (rather than just zero or per- 
fect bonding), various levels of shear resistance, and rational methods of 
incorporating the weight and strength of the glue joints and exterior cover. 
(Note that any SEA changes due to glue joints or freezing of the crushable 
material could be incorporated by changing fig. 5 and the corresponding 
equations defining the material.) 

Two alternate models have been briefly investigated; a full investigation 
would require changes in the machine programs or the design charts. In one, 
called the "hemisphere" model, the payload and the entire upper hemisphere are 
able to penetrate as a unit, regardless of bonding, because of low resistance 
to cross-grain crushing in the equatorial plane. In the other, called the 
"shear-plug" model, the payload and the cylinder (with rounded ends) directly 
above it are able to penetrate as a unit, again regardless of bonding, because 
of low shear resistance over the cylinder walls. 

The investigation of the hemisphere model indicates that the payload and 
upper hemisphere would have started to penetrate below the maximum g load- 
ings for almost all cases in table 1, and would have done so more readily than 
the isolated payload for the detailed and simplified models, provided that the 
cross -grain crushing strength is 8 percent of the end grain value. Hemisphere 
penetration is calculated even at 18 percent (ref. 3 for balsa) for several 
cases. At 8 and 14 percent, however, it is also calculated for the high-speed 
(374 ft/sec) experimental configuration at a g loading of 3000 to 4000 (ref. 
10); yet there was no evidence of hemisphere penetration (i.e., of exterior 
cover wrinkling at the equator). Thus, 18 percent may be the most nearly cor- 
rect figure, particularly when the glue joints are considered. In any event, 
the hemisphere model is of interest. 

For the shear-plug model, on the other hand, there is no possibility of 
shear-plug penetration at the g loadings under consideration according to 
calculations based on a measured shear strength. This model would be of 
interest only for much higher g loadings (say, greater than 3000) and/or 
much lower shear strengths (less than 145 psi) . 


22 


I 



Finally, a question remains as to the effect of payload radius Rp when 
densities (including payload packaging density) and stresses remain the same, 
that is, a question as to the effect of scaling. Except for the trivial grav- 
ity force, the equations of motion indicate (for a nondimensionalization dif- 
ferent from that of appendix A) that geometric scaling should apply, with U Q 
and SEA unaffected, with n p max ~ 1/Rp, R ~ Rp> and w co ~ R p 1 * 3 * These scaling 
conditions, then, can be used to extend the present results. They can also be 
used to check numerical results whenever two cases scale each other. In the 
present examples, only cases 3 and 9 should roughly scale each other, as 
indeed they do. 


CONCLUDING REMARKS 


The governing equations for the landing impact of a rigid payload 
protected by a crushable casing, including the possibility of penetration by 
the payload into the casing, have been developed for a general vertically 
symmetrical landing. The general equations have been specialized for zero 
shear resistance, a constant compacting strain, and uniform density of the 
crushable casing. They have also been specialized for a spherical payload and 
casing, with the latter having its highest crushing stress in the radial direc- 
tion. For the spherical configuration, two approximate analytical models have 
been defined: (1) a detailed model with no additional assumptions but requir- 

ing an automatic computer program; and (2) a simplified model utilizing either 
the computer program or design charts but requiring the assumptions of infi- 
nitely thin sheets of crushed material, constant mass in the equations of 
motion, zero acceleration due to gravity, and a constant resistance to 
payload penetration. When specialized to prevent penetration, the simplified 
model (or slight variations thereof) has been widely employed in previous 
work; and the specialization of the detailed model for a uniform and isotropic 
crushing stress without penetration has been shown to have the same basic equa- 
tion as an earlier model. Results for the two models are in reasonable agree- 
ment with two previous measurements, having impact velocities of 220 and 
374 ft/sec. 

Thirteen examples have been presented for each of the two analytical 
models. The examples are for an impact velocity of 300 ft/sec, a maximum 
permissible g loading of 2000, payload weights of 100 and 450 lb, payload 
radii ranging from 0.6 to 1.6 ft, payload packaging densities ranging from 
26.23 to 110.5 lb/ft 3 , a choice of zero or "perfect" bonding between the pay- 
load and the crushable material, and a choice of a selected balsa-like or 
honeycomb- like class of crushable material. Overall radii resulting from the 
designs range from 1.81 ft to about 3.45 ft, and the resulting crushable 
material weights vary from a little over 124 lb to almost 763 lb. The fol- 
lowing conclusions are drawn from the examples: 

1. The simplified model has the lower crushable material weight of the 

two models, the greatest difference being approximately 15 percent 

for an example without payload penetration. 
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2. The greatest difference between the two models is a 29 percent 
discrepancy in the maximum crushing stress determined by a design 
example without penetration for the honeycomb- like class of material. 

3. In the absence of payload penetration, the crushable material weight 
and the overall radius increase with increasing payload radius, that 
is, with decreasing payload packaging density. 

4. For the honeycomb- like class of material and the 100- lb payload (with 
the maximum g loading held at 2000 to minimize weight) payload 
penetration does not exist for the detailed model. For the simpli- 
fied model, however, penetration slightly reduces both the crushable 
material weight and the overall radius at the lowest payload radius, 
0.6 ft, the only radius among those tried for which penetration 
exists with honeycomb. This advantage of penetration means that 
bonding should be avoided with the present honeycomb -like material, 
for which the specific energy absorption increases with crushing 
stress (a variation more common even for balsa than that for the 
present balsa-like class). 

5. For the balsa-like class of material (with a payload weight of 
450 lb), penetration drastically increases the crushable material 
weight at the lowest payload radius (1.0 ft) and slightly increases 

it at the next lowest radius (1.2 ft). This combines with nonpenetra- 
tion results to give a minimum-weight radius for an unbonded payload 
between 1.2 and 1.6 ft and a corresponding minimum-weight payload 
packaging density between 26.33 and 62.15 lb/ft 3 . 

6. An impressive benefit of payload penetration is a decrease in 
crushable material weight by a factor greater than 4, which occurs 
for the 100-lb payload with a radius of 0.7 ft when the (selected) 
honeycomb -like class of material is required without penetration (by 
the g loading ceiling of 2000), but when the more efficient balsa- 
like class is feasible with penetration. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, Calif., 94035, Jan. 12, 1970 
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APPENDIX A 


DEVELOPMENT OF BASIC EQUATIONS 
GENERAL GOVERNING EQUATIONS 


Figure 2 shows a fairly general landing geometry at the start of impact 
and at two later instants of time during impact. The only shape requirement 
is symmetry about two mutually perpendicular vertical planes, and this 
requirement is made for compatibility with the assumption of pure vertical 
translation. The crushing strength and density of the crushable material are 
also required to be symmetrical but otherwise generally variable. The land- 
ing surface is assumed to be perfectly flat and rigid. The payload is 
assumed to be perfectly rigid and is considered unbonded to the crushable 
material in determining the start of relative motion (if any) , herein 
referred to as M payload penetration.” 

The analysis begins with the phase of impact shown in figure 2(b). In 
this phase, the crushable material has begun to crush against the landing 
surface, but payload penetration has not begun. Hence, 

q p = q (ai) 

where qp and q are the absolute displacements of the payload and the 
uncrushed crushable material, respectively, as shown in figure 2(b). The 
equation of motion for the constant payload mass mp 0 is written as: 

m q = m g - / / a dA 

po 4 po^ J J vpo (A2) 

A po 

where q is the second time derivative of q, g the local acceleration due 
to gravity, a v p G the vertical component of normal stress on mp Q , and Ap 0 
the area over which the stress acts. 1 

The equation of motion is now written for the variable mass m co - m cl 
of uncrushed crushable material in figure 2(b) under the assumption that each 
particle in m co - m c i is moving at the same vertical velocity q at a 
given instant. The equation is 

x Note that the vertical component of shear stress on the payload could 
have been added to o vp 0 in equation (A2) and in the later equations con- 
taining crypo. There are, however, no aypo terms in the major governing 
equations except for the equation that determines the existence of payload 
penetration. When penetration is about to start, the most important stresses 
are aypo values that are almost large enough to crush the material (in its 
strongest direction if it has one and has been deployed to utilize the fact) . 
Shear stresses are negligible by comparison and are often incorporated in the 
tests to determine the crushing strength. 
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( m co - ra ci)q = C m 


Cl- 


CO 


m )g + 


f f\ po dA - / / 


a dA 
v 


(A3) 


po 


xi 


where a y is the vertical component of the "mostly static" and "mostly 
normal" stress capable of deforming the crushable material plastically oyer 
the area A C1 . The term "mostly static" implies the assumptions that the 
damping force is incorporated and that the coupling between vertical and 
horizontal velocity, or dynamic buckling effect., is also incorporated. The 
term "mostly normal" implies the incorporation of a small shear stress. (See 
the discussion of figs. 1(b) and 1(c) under "Properties of Typical Crushable 
Material.") The area A C1 is the intersection between m co - m C i and the 
mass of crushed material m cl in figure 2(b). The crushed material is 
assumed to have been compacted sufficiently to deliver the required crushing 
stress a v ; and successive layers of material are transferred from m co - m cl 
to m cl at the boundary A cl as they undergo plastic deformation. (Such a 
process is observed experimentally in the crushing of plastic foam, balsa 
parallel to its grain, and honeycomb parallel to its axis.) 

The assumption of a uniform (though time variable) vertical velocity q 
throughout m c0 - m cl has two implications besides the crushing by layers 
just described. First, the elastic stress waves that establish a uniform q 
must obviously travel at speeds far greater than q (i.e., q must be a low 
subsonic value) . Second, the deformation waves associated with the elastic 
stress waves must be small enough not to affect the magnitude of q. 

Equations of motion could also be written for m C 2 in figure 2 and for 
the elemental layer of mass being transferred to m c i , These equations 
would determine dynamic stresses in the sense of variable mass, however, and 
are not needed since a v is assumed to be a measured quantity. 


Equations (A2) and (A3) can be used to determine whether payload 
penetration occurs and, if so, at what displacement. If both equations are 
solved for q - g and the results equated and rearranged, the result is: 


//■ 


dA 

vpo m. 


^o 


^po 


+ m 
po co 
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jj 


a dA 
v 


(A4) 


v ci 


where nip 0 + m co - m cl is the total time-dependent mass at velocity q in 
the absence of payload penetration. Equation (A4) is the only major governing 
equation containing ayp 0 ; the neglect or incorporation of shear stresses is 
justified in footnote 1 of this appendix. 


If equation (A4) is considered formally solved for aypo, payload 
penetration cannot begin until those stresses become large enough to cause 
plastic failure, that is, become a a v distribution. This will occur only 
if the geometry, masses, and impact conditions are such that the displacement 
q becomes large enough to bring k c i up to the necessary size. It can be 
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seen in equation (A4) that the following quantities tend to prevent or post- 
pone payload penetration: a large payload bearing area Ap 0 , a small payload 

mass nipo, a large total mass mp 0 + m co - m c i, and a small bearing area of 
compacted material A cl . 

Equations (A- 2 ) and (A3) can be added to give the major governing 
equation: 


(m 


po 


+ m 


co 


- m 


, ) q = (m + m 
v po 


co 


- m 


^ ■//' 


dA 


(A5) 


v cl 


This is the simplest equation of motion to apply during the impact period 
prior to payload penetration. It does not contain crypo (see footnote 1 of 
this appendix) . 


If equation (A4) shows that payload penetration has begun, the geometry 
of figure 2(c) applies. The constant payload mass nip 0 is then positioned 
by the coordinate qp so that equation (A2) becomes: 

mq=mg-//a dA 
po M p po & J J u ypo (A6) 

A po 

The mass ^pi in figure 2(c) consists of compacted crushable material. 
The velocity qp of mp G is assumed to apply uniformly throughout . 

Hence, the equation of motion for m-pi can be written: 

m, 5 = m g + / / a dA-//a 1 dA 

P 1 T PI 7 J V P° J J v (A7) 

A A . 

po pi 

where Ap 1 is the intersection between m-pi and the variable mass 
m co " m cl ” m pl °f uncrushed crushable material, and o v l differs from the 
"mostly static" crushing vertical component a v acting over Ap X only 
because the elemental layer of mass dmp^ is transferred to mpi . 

The equation of motion for dmp X becomes, with the definitions just 
given: 


or 




dA 




(A8) 
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where mp X = dmp X /dt, where the higher-order quantity gdmp X has been 
neglected, and where qp and q represent the final and initial velocities 
of dmp X . 

Equations (A6) , (A7) , and (A8) can be added to give 

* V^p * % - * ( v * V ,g -f A dA (A9) 

Apl 

A second major equation of motion defines the forces on the variable mass 
m co " m ci ” ropi °f uncrushed crushable material in figure 2(c). Under the 
assumption that each particle in m co - m cX 
vertical velocity q, this equation is 


(”co - m cl - V 1 ' 4 * (m co - m cl - V g 


where Ap X is defined as in equations (A7) 
analogously to that in equation (A3). 

Equations (A9) and (A10) define the motion during payload penetration 
since all quantities can be defined in terms of q and qp and their deriva- 
tives. The two major governing equations are coupled in the most general 
case by the qp - q term, the dependence of mp X on qp - q, and the 
dependence of nip X and of o v over Ap X on qp - q. Equations (A9) and 
(A10) do not contain a-ypp. 


mp X is moving at the same 


*/ /° » " ( A dA 

Api A Cl 


(A10) 


and (A8) and A cX is defined 


GOVERNING EQUATIONS FOR ZERO SHEAR DEFORMATION, 
UNIFORM COMPACTING STRAIN, AND UNIFORM 
MATERIAL DENSITY 


Although figure 2 does not show any shear deformation, the equations 
derived so far would apply even if such deformation were present. The speci- 
fication of a v , however, requires knowledge of the location of the surfaces 
A cx and Ap X , and shear deformation is ruled out at this point to retain the 
simple surfaces implied by figure 2. Then A cX and Api are determined by the 
height variables h cX and hp X , where h cX is the local height of m cX in 

figures 2(b) and 2(c) and where hp x is the local height of m pl in fig- 
ure 2(c). If shear deformation had been considered, there would be a trail- 
ing of material to make mp X wider than mp Q and a possible lifting off 
the ground of the edges of m cX . 
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The actual lengths of hp^ and h c ^ are evaluated in terms of a readily 
measurable variable e. This variable is defined as the compressive compact- 
ing strain of the crushable material ("compacting strain" is the strain at 
which the crushing stress rises abruptly from a relatively constant value for 
a test specimen of uniform cross section) . The effect of the direction of 
maximum strength is ignored under the assumption that the material will com- 
press to the same compacted strain along any axis. 


The use of e is made feasible by the assumption that the energy 
absorbing process consists of the vertical crushing of separate vertical rods 
of material. Then h c in figures 2(b) and 2(c) is the total shortening 
deformation of a rod having an initial length of h c + h C i, and qp - q is 
the corresponding deformation of a rod having an initial length of 
qp - q + hp! , where hpi is the local height of mpi in figure 2(c). 

Since the variable e dl is (by the definition of strain) the total shorten- 
ing deformation of each successive rod element d& to be crushed, it is seen 


that 


A 


q D - q 


-f 


q p -q+hpi 


dfe 


(All) 


h c 



cl 


e dl 


J 


Under the assumption of uniform compacting strain (e = = constant), 

the equations just derived can easily be solved for h-pi and h c i 


hpi 



h 


cl 





J 


(A12) 


The effects of equations (A12) for the constant can be seen in figure 3, 

which is otherwise identical to figure 2. In figures 3(b) and 3(c), the mass 
m c i has become a foreshortened image of the volume V c , which would lie 
beneath the landing surface if there had been no crushing; and the height 
hpi of the mass nipx has become constant in figure 3(c) . 


The next limitation to be imposed in the analysis is that of uniform 
density of the crushable material prior to crushing. Then mpi and m c i in 
equations (A4) , (AS), (A9) , and (A10) become: 
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m = p . (V + V ) 
pi ck pi p' 


m 


Cl - p ck< V cl * V 


J 


where V pl is the volume of m pl , V p the volume swept out by m po during 
payload penetration, and V C i the volume of m cl . When V pl is related to 
V p and Vci to V c according to equations (A12) , the masses m pi and m C i 
become 


1 - £, 


m = p , IV 

pi ck\ p e 


+ V T 


m cl p ckfc e. 


+ V 


P ,v 

ck p 


p ck V c 


J 


(A13) 


The relations between the volumes V p and V c and the displacements q 
and q p can be determined by the following volume formulas: 

'\ 


V = A . (q - q) 
p poh ^ H p n J 


-r 


A , dh 
sh c 


CA14) 


J 


where Ap 0 ^ is the horizontal planar projection of Ap 0 and Apj, and A s h is 
a horizontal cross section in V c . 


Another useful set of relations involves the stresses. The stress 
integrals in the major governing equations (A4) 9 (AS), (A9) , and (A10) can be 
written 






a dA 
po 


po 


poh 


/ A dA f f ° dA 


Cl 


v clh 




a dA 


poh 




J 


(A15) 
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where 0p O is the normal stress on the payload prior to penetration, o is 
the normal crushing stress on any surface, and A^h is the horizontal 
planar projection of A C i . For integration of the right-hand sides of equa- 
tions (A15) , the normal stresses must be evaluated at the intersections 
between the curved surfaces (Ap 0 , A cl , Apj) and vertical lines through the 
centroids of the elements dA in the horizontal planar surfaces (Ap 0 h> A C ih, 
A poh) * 

A final useful relation is the following change of variable: 


With equation (A16) 



- _ dq _ dU 1 1 dU 1 d(U 2 ) 
q = dt = dt = U Hq = 2 — dq" 


(A16) 


(A1 7) 


Equations (A13) for uniform density and uniform compacting strain, 
together with the volume equations (A14) and the stress equations (A15) , are 
now introduced into the basic governing equations. In addition, one of the 
governing equations is reduced to first order by equation (A17) . Thus equa- 
tion (A4) , which determines the start of payload penetration (if any), 
becomes : 



0 

Equation (A5) , which defines the impact prior to penetration, becomes: 


1_ d(U 2 ) 
2 dq 



Finally, equations (A9) and (A10) 


for payload penetration become: 


(A18) 


(A19) 
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q p = g - 


(P^A^/evKi, - q) 2 fa dA 

poh 


*\ 


ck poh' k J L4 p - A 


q = g 


m + (p ,A , /e, ) (q - q) 
po ck poh k n p n 

f fe dA -/ fo dA 
^poh 




CA20) 


cih 


>/ 


m co - ^ p ck /e k^ / A sh dh c - (p C k A p 0 h /e k ) (q P ~ q) 


J 


Equations (A18) through (A20) define the problem for the landing geometry of 
figure 3, providing the stress and volume integrals can be evaluated. 


EVALUATION OF VOLUME AND STRESS INTEGRALS 
FOR SPHERICAL GEOMETRY 


The evaluation of the volume and stress integrals is facilitated by 
restriction of the landing configuration to a simple shape. The sphere is 
particularly useful for landing packages and crushable coverings that must 
absorb energy from impacts in any direction. Hence, the geometry is now 
specialized to concentric spheres. 

The spherical landing package is shown in figure 4(a) at the start of 
impact, in figure 4(b) during impact but prior to payload penetration, and in 
figure 4(c) after penetration. All angles and radial lengths are shown in a 
plane of maximum value. 

The volume integral in equations (A18) through (A20) can be evaluated 
immediately on the basis of sketch (a), which is a simplification of figure 



Sketch (a) 
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4(a) or 4(b). Thus 


or 



(R - q + h c ) 2 ]dh c 



j q 2 (3R - q) 


(A21) 


It is also noted, for certain volume terms and integration areas in equa- 
tions (A18) and (A20) , that 


Apoh - v Rp 2 


(A22) 


as can be seen from figure 4 and the definition of Ap 0 h* 

One of the stress integrals can be evaluated immediately under the 
temporary assumption that expo = crpok = constant. Then, with equation 
(A22) , 


/ / a dA = a , A , = ttR 2 a , 
J po pok poh p pok 

A poh 


(A23) 


The assumption that a p0 = a po k ma Y be valid only when a po k approaches 
the normal crushing stress a, that is, when payload penetration is 
approached. Fortunately, this is the only region of interest for 
equations (A18) and (A23) . 

Prior to evaluation of the stress integrals for a, another limitation 
is applied to the analysis, namely, 

a = a Q cos a ; a < 90° (A24) 

where a is the angle between the local normal to the stressed area and a 
radial line from the center of the spherical system (as if undeformed) through 
the point of stress application. (The rationale for this limitation is given 
immediately prior to eq. (2) .) 

The next step is to determine a for m c i in figure 4(b) (with 
applicability after payload penetration as well) . Any point on the surface 
of revolution shown by the dotted portion of the circle is defined, with the 
aid of sketch (a) , by 

r 2 + (R - q + h c ) 2 = R 2 
From the second of equations (A12) 
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(A25) 


r 2 + I R - q + 




R 2 


l - e . n cl 
\ 1 e k 

Hence, in the vertical great-circle plane, the surface slope of m cl is 


dh 


cl 


dr 


__ u - •k ,r 

E k( R ■ q * [ (t k /(i ■ E k j ] hci ! 


(1 - e k )r 


e, / k 2 - r 2 
k 


Thus, with (J> shown in figure 4(b) and from the greater detail in sketch (b) , 



Sketch (b) 
dh 


a = <b - arc tan 


ci 


dr 


+ arc tan 


a - £ k )r 


e^/R 2 - r 2 


(A26) 


With equation (A26) introduced into equation (A24) , the normal stress on 
the surface A cl (at height h c i in fig. 4(b)) becomes: 


o=o cos 
0 


Thus , 


<f> + arc tan 


U - £ k )r 


e^/fo 2 - r 2 


a (1 " £ k )r 

— = cos 4> cos arc tan 

o 


A 2 - 


- sin d> sin arc tan 


(1 - e k )r 

e, / r 2 - r 2 
k 


but, in figure 4(b) and sketch (b) , 


R - q - h 


cos <j> = 


Cl 


/r 2 + (R - q - h ) 2 
^ Cl 


sin ip = 


/r 2 + (R - q - h ci ) 
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and h c i can be determined from equation (A25) . Then, after some 
manipulation, together with 


we obtain 



VR 2 -(R-q) 2 


a 2irr dr 


/r 2 - (R-q ) 2 


(/“‘"I ( 

n ' 


[ (1 - e-O 2 * 1 


2tt[(R - q)/R 2 - r 2 - (1 - e k )R 2 ]r dr 
e k 2 (R 2 - r 2 ) ] Jr 2 + (l/e k 2 )[R - q - (1 - /R 2 - r 2 ] 


With the substitution b 


/ R 2 - r 2 , this becomes 



a dA 


A 


clh 




[(R - q)b - (1 - £ k )R 2 ]b db 

(1 - e k ) 2 ]b 2 ♦ (1 - E k ) 2 R 2 ]|R 2 - b 2 + (l/e k 2 )[R - q - d 



(A27) 


and equation (A27) is the stress integral for the upper surface A cl of 
m cl in figure 4(b) or 4(c). The integral contains no singularities and 
hence can readily be evaluated numerically. 


The angle a is now determined for mpj in figure 4(c). Thus, with 

sketch (c) for detail, 
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or, with the first of equations (A12) , 


a = arc sin — arc tan 

K-n 


d/e k ) Cq p - q) + A p 2 _ r p 2 


(A28) 


When equation (A28) is introduced into equation (A24) , the normal 
crushing stress a is determined. When a is put into J^ Q ^f 0 dA with 


Rr 


• r Ap„J 0dA *( ^P^P, 


the result is 


// 

A u 
poh 


a dA = c 


2 tt I” Cq - q)/e, )/r 2 - r 2 + R 2 lr dr 
L H P 4 ^ P P P J P 1 


o P 


R„/ 2[Cq„ - q)/ £ k ]v/ R p 2 - rp 2 + [(q T , - q)/e v ] 2 + R„ 2 


With the substitution b p = /R p 2 - r p 2 , this becomes 


// 


2t ra 


a dA 


| [ (q - q)/e,]b + R 2 )b db 

I t h k J P P I p i 


A 


poh 


/ 2 [(q - q)/e k ]b + [ (q - q )/^] 2 + 


(A29) 


Equation (A29) is the stress integral for the lower surface A pl of m pl in 
figure 4(c). It can be integrated directly to give 


If 

A u 
poh 


2tt<7 


a dA 


ISR^ 


q P - q 


p 


-(V) •( 


4 ,2 

q p ' q 


R p 2 ♦ 3R p ^ 




&r) - ■ 


3Rn 


q P - q^ 


‘P 


(A30) 


Since equation (A30) is indeterminant at [ (q p - q)/e] = 0, that is, at 
the start of payload penetration, the best way to evaluate / Vh / a dA is to 

integrate equation (A29) numerically. At [Cqp - q)/e] = 0, equation (A29) 
becomes 
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(A31) 



Comparison of equation (A31) with equations (A23) and (A24) serves as a 
partial check on the development of equations (A29) through (A31) since 
a = 0° (or a = a Q ) on the payload at the start of penetration. 


Dimensionless Governing Equations for Spherical Geometry, 
With Termination Conditions 


Governing equations that have been derived in this appendix are now 
made dimensionless in a form convenient for numerical solution of the impact 
problem with spherical geometry. Prior to payload penetration, equation 
(A19) applies and becomes, with equations (A21) and (A27) , 


n P _ 1 d(w 2 ) g „ Fcl(z 
n md 2 dz ' n md§e R m(z) 


(A32) 


where g e is the acceleration due to gravity on earth, npg e and n mc jg e are 
the actual and maximum design payload decelerations, respectively, and 


Id „ 2 ^ 

Ve ' - 2 d? U 


w = 


A 


n mdSe R 


Z 5 R 


ttR 2 o 


K r = 

K m. 


po^idSe 




J 


(A33) 


The dimensionless total mass m(z) in figure 4(b) is defined by 


m(z) = 


m po + m CO - m cl 


m 


m. 


= 1 + 


co 


po 


m po ^ e k 


1_ / p ck\ 
e k \ p pR/ 


z 2 (3 - z) 


(A34) 


with 


m 


l po 


pR (4/3) ttR 3 


(A35) 
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where F cl (z) is the dimensionless crushing force on m(z) defined by 


Foltz) 



2 


/ 


1 



(1 - e k ) 2 ]s 2 + (1 


[(1 - z)s - 1 + e k ]s ds 
- e k ) 2 }(l - s 2 + (l/e k 2 )[l-z- (1 



(A36) 


with s = b/R. Under the initial condition w(z = 0) = w Q , equation (A32) 
integrates immediately to 


w 


2 



w r 


+ 


+ 



FciCj)] 

m(j) 


= w G 2 + 


2gz 

n md®e 


2Kr 



z F c i (j) 
m(j) 


d J 


(A37) 

Equations (A32) and (A36) apply until z reaches a value z s , at which 
payload penetration (if any) starts. The start of penetration is defined by 
equation (A18) with equations (A23) , (A34) , and (A36) . Thus, 

°pok _ /r \ 2 F dCz) 

a o " \Rp/ m(z) CA38) 

and penetration starts (z = z s ) if and when Opok reaches a G (providing the 
payload is not attached to the crushable material). 

At the start of penetration, it becomes convenient to introduce three 
new variables and one constant as follows: 
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"A 


q R 

y E R^ = Z E 


x = t 


' n md®e 


*P 


e = 


q p - q 

e k R p 




_ V°o _ /V 

P n ’po n md2e R \ r ) 

The initial conditions for penetration are 


J 




r o', — — R 

Y(x = 0) = y s = z s 


^dy\ = /dy N 

d x L " H ‘ " s 


w, 


flT 

*P 




e(x = 0) = e,. = 0 


( 


SfL- si- 


J 


(A 39) 


(A40) 


where w s is the value of w determined by equation (A37) at z = z s . 

During penetration, equations (A20) apply, with equations (A21) , (A22) , 
(A29) , (A36) , and (A39) . Thus 


n 


P 

n md 


s d 2 e + 

£ k “T + 

dx 


d 2 y 

dx 2 


(iTRp 3 Pck £ k/ m po^ (de/dx) 2 + KpFp Q (e) "A 


n md®e 


m p en(e) 


> (A41) 


d 2 y = g 

dx 2 n md®e 


K R R ciCy) - K p F po (e) 
m cr (y,e) 
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where F cl (y) is determined by equation 036) with (R p /R)y substituted for 
z; m pen (e) is the dimensionless total penetrating mass in figure 4(c) defined 
by 


m pen ^ 


"Vo 


+ m, 


l Pl 


m- 


= 1+4 


po 



(A42) 


with 


m po 

P P = T 

(4/ 3) nRp 

m cr(y> e ) is the dimensionless crushable mass in figure 4(c) defined by 


(A43) 


m (y,e) 
cr w * 


m co m ci m pl m co 
^o m po 



(A44) 


and Fp 0 (e) is the dimensionless crushing force on mp en (yp,y) defined by 


F 

po 


(e) 


with Sp = bp/Rp. 


" R p 0 o a. 


// 


0 dA = 


poh 



(es p + l)s p dSp 
y^eSp + e 2 + 1 


(A45) 


Equations (A32) , (A37) , (A40) , and (A41) are the required governing 
equations and initial conditions. With rebound excluded and with no payload 
penetration, the impact is terminated (z = Zmax) when w = 0. Thus 


w ( z max) = 0 


(A46) 


With rebound excluded but payload penetration present, the impact is termi- 
nated (y = y ma x y + E k e = (y + e ke) max ) when dy/dx = 0 and 
dy/dx + ek(de/dx) = 0. Since e increases monotonically with time for the 
present examples, the termination condition becomes 



(A47) 


As the payload begins to penetrate the crushable material, equations 
(A41) are integrated simultaneously. This continues until the first of 
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equations (A47) is satisfied, thereby providing initial conditions for the next 
phase of the problem. In this next phase, the second of equations (A41) can 
be ignored and the first equation solved alone. The uncoupled solution of 
the first of equations (A41) continues until the second of equations (A47) is 
satisfied. 


Simplifications for Constant Mass, Infinitely Thin Crushed Material, 
and/or Constant Penetration Resistance 

The governing equations are now written with constant mass (p c k = 0) . 
Equations (A32) and (A34) are combined to give 


1_ d(w 2 ) = 
2 dz 


- % 


FciCz) 


(A48) 


n md - - "md & e 1 + (m^/m^) 

where F c i (z) is given by equation (A36) . Equation (A37) becomes 

Fci(j) 


w 2 = w 0 Z + 21 


2 f 

J 2 d 3 J o J n md g e 
0 0 


- K 


R 1 + (m co /m po ) 


dj 


= "o 2 * ^ 


2Kr 


n md^e 

and equation (/\38) gives 


1 + (m /m 
v co' po 




F c i(3)dj 


(A49) 


J pok 


,R_ 

Rr 


F c iU) 

1 + (n, co /m po T 


(A50) 


After payload penetration starts, equations (A41) apply, and p c ^ = 0 gives 


n r 


n 


= e. 


md 


d^£ + £X. = .. g . . k_F (e) 

dx 2 dx 2 We > P° C 


d y g_ 


^d^e 

K R F c i(y) - KpFpo(e) 


dx" 




“co^po 


(A51) 


where F c i (y) is given by equation (A36) with y(Rp/R) substituted for z and 
where Fp 0 (e) is given by equation (A45) . 

Regardless of whether p c k = 0 or ^ 0, equation (A36) is greatly 
simplified by the specialization to = 1 (infinitely thin crushed material) . 
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Thus, for efc = 1, equation 036 ) can be integrated; and z = y(Rp/R) can be 
substituted as follows: 


F cl (z) = 2z (1 - z) ] 



052 ) 


Regardless of whether = 0 or / 0 and whether - 1 or / 1, equation 

(A45) is greatly simplified by the assumption that the payload stress and 
force maintain their initial penetration values throughout the penetration 
stroke, that is, the assumption that Fp 0 (e) can be represented by Fp O (0), 
which is unity. Thus, the assumption has the form 

FpoCe) = F po (0) = 1 053) 

Equations (A52) and (A53) can be used individually or together in the govern- 
ing equations for p c k / 0 or for = 0 whenever the simplifying 

assumptions are appropriate. 

If the assumption of equation (A53) is made (constant penetration 
resistance) and if = 0 (constant mass), equations (A41) can be decoupled 

by a change of variable, regardless of whether = 1 or / 1, as can equa- 
tions (A51) , which are specialized for p c k = 0. The governing equations 
(A48) through (A51) will be rewritten, however, for = 1, that is, equa- 
tions (A52) , as well as the assumption of equation (A53) . The change of 
variables is 


y 


q_ = r_ 

Rp hp 




v 


dq/dt 

/n mdg e R p 


U 

/n mdge R p 



qp 

y p E E e k e + y = e + y 




_ dq P /dt _ de dy = de dy 

V P ~ f — dx + dx dx dx + dx 

*Wp J 


(A54) 


With the first and second of equations (A54) and the first of equations 
(A52) , as well as equation (A39) , equation (A48) becomes 


Hp _ 1 d(v 2 ) = g 2K y[(R/R P } " y] 
"md " 2 dy n md g e P 1 + 


(A55) 
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Equation (A55) integrates immediately to 

.y 

i a ( ■>r2.'\ 

d? = v 0 


V 2 = Vo 2 + 2 / i d£ = v „2 + __?JQL_ . f V ^ R/ V (2/3)y] 


J 1 d(v 2 ) 

7 2 d£ 

* / 0 


n md®e 


1 + (m /m ) 
v co per 


(A56) 


The same substitution in equation (A50) gives 

a p ° k / R \ 2 2 (Rp/R)y [l - (R p /R)y] 


Rr 


1 + (m co /mp 0 ) 


(A5 7) 


and equations (A55) and (A56) apply until payload penetration starts at 
CT pok/°o = 1 in equation (A57), which then becomes 


y s 2 - 





0 


CA58) 


where y s is the dimensionless displacement at the start of penetration. 

After penetration begins, equations (A51) apply. They become, with 
equations (A52) , (A53) , and (A54) , 


np_ i d(v p ^ „ 

n md 2 dy P n md 8 e P 


1^ d (v 2 ) 
2 dy 


g 

n md^e 


K p j2y[(R/R p ) - y] - 1 


m /m 
co' po 






(A59) 


Equations (A59) are uncoupled. The first equation integrates immediately to 


V P 


2 


= V s ‘ 



(A60) 


and the second to 


v 2 = v s 2 + 2 


n ,g m /m 
md & e co po 


(y + y s ) + I (y2 + yys + y s 2 ) 1 o - y s ) 

(A61) 


where v s is the dimensionless velocity corresponding to y s . 


Equations (A55) through (A61) are the governing equations for infinitely 
thin crushed material = 1) and constant mass (p^ = 0), with the assump- 
tion of constant penetration resistance (Fp 0 (e) = 1) required after payload 
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penetration but not before. In the absence of penetration, equation (A56) 
can be used to determine y max by setting v = 0. Thus y max is the solu- 
tion of the cubic equation 


2K P y ,L[wV - < 2 ^r ma x] 

1 + (m co /m p o) 


2gy. 


max 


n md g e 


- v, 


2 _ 


= 0 


(A62) 


The dimensionless velocity v s at which penetration starts (if any) can be 
determined by writing equation (A56) as 


v 2 = v 2 + 
s ° n 


2gy s 2K p y s 2 [ (R/Rp) - (2/3)y s ] 


'md®e 


1 + 


(m /m ) 
v co po' 


(A63) 


where y s is the solution of equation (A58) . Then the maximum payload 
dimensionless displacement y p with penetration is determined by equa- 
tions (A60) and (A63) with Vp m = X 0 as 


Vo 2 - 


y* 


max 


2K p y s (l - jy s [(R/R p ) - (2/3)y s ]/[l + (m co /m po )]} ) 

2[K p - (g/Hmdge^ ~ 


(A64) 


A similar procedure with equations (A61) and (A63) and with v = 0 gives 
ymax as the solution of the cubic equation. 


2K p 

m /m 
co po 


2 r 2 

(y + 

3 max 


y y 

'max s 


^s ^ Rp ^max + ^s^ 


+ 1 


cy, 
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y s ) 


(y - y ) 
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2gy. 


max 


n md^e 


+ v r 


2Kpy s 2 [(R/R p ) - (2/3)y s ] 


1 + 


(m /m ) 
co po' 


= 0 


(A65) 

It should be noted that yp max is generally sufficient for design problems 
involving payload penetration and that equation (A65) for ymax is included 
only for completeness and possible checking purposes. 
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APPENDIX B 


SAMPLE CALCULATION FOR DESIGN BY SIMPLIFIED MODEL 
WITHOUT PAYLOAD PENETRATION 

Under the assumptions that Pck = 0, e k = 1, and g = 0 (simplified 
model), a calculation is performed in this appendix to yield a design of a 
crushable casing for the following conditions: 

U Q = 300 ft/sec 

ed = 0 . 7 

e m ~ 0*8 

n D = 2000 
Umax 

w po = mpo^e = 100 lb 
Rp = 0.6 ft 

Payload penetration prevented by bonding 
if necessary 

The calculation corresponds to case 1 in table 1, as indicated by the last 
four of the above conditions. The calculation follows: 


Uo 2 


(300) : 


( n p max g e )Rp (2000) (32.17) (0.6) 
From figure 8 for = 0.7 


= 2.331 


= 0.286 


or 


R = 


0.6 

0.286 


= 2.097 ft 


From figure 9 for = 0.7 with Rp/R = 0.286, 

J ma = 2.74 


check : 


Uo 2 


(n Pmax^ R 


= (2.331) (0,286) = 0.6665 


From figure 6 


"max 


hnax 

R 


= 0.500 
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or 


q = 0.500(2.097) = 1.049 ft 
M max 

From figure 7 for = 0.7 

j ma =2.74 (check is perfect to 3 places; therefore use 2.74) 
From figure 12 with Rp/R = 0.286, 

Njna 55 55.5 

e m^ma 1=1 0.8(2.74) = 2.192 


ge £ d 


32.17C0.7) 
(300) 2 


2.503xl0 -4 ft" 1 


JmaWpo 
N mo TrR p 3 (144) 


(2.74) (100) 

(55 . 5) tt (0 .216) (144) 


0.5053X10" 1 psi/ft 


From equation (27) with a Q in psi 


SEA = — 2 - 192 — ft- lb/lb 

2.503xl0" 4 - (O.SOSSxIO-Vao) 


Now decide which material to use in figure 5 (trial eliminates the balsa- like 
class, which is too strong for low payload weight). From figure 5 for 
honeycomb- like material, 

try SEA = 10,400 ft-lb/lb at a Q = 1,000 psi 

SEA - ■ 10 - 960 ft - lb/lb 
try SEA = 10,870 ft-lb/lb at a Q = 1,100 psi 

SEA ° 2.5of-*0 P 45 9 - 10 - 72 ° ft - lb/lb 
try SEA = 10,800 ft-lb/lb at a Q = 1,090 psi 

SEA = 2 . 1 4^0 63 = 10 ’ 750 ft ' lb/lb 
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try SEA = 10,770 ft-lb/lb at a Q = 1,080 psi 
2.192xl0 4 


SEA = 


= 10,770 ft-lb/lb 


therefore 


2.503 - 0.468 
a Q = 1,080 psi 
SEA = 10,770 ft-lb/lb 

From equation (8), with o Q in psi and SEA in ft-lb/lb 


P £ 
cm s e 


£ m( 14 4 Pq) _ (Q.8) (144) (1,080) _ 


SEA 


10,770 


= 11.55 lb/ft 3 


From equation (17) 
m 


co 


m po + m co 


J ma U o 2e m 

g e (SEA)e d 


(2.74) (300) 2 (0 . 8) 
(32.17) (10V7'70) (0 . 7) 


= 0.8135 


"Vo 


m, 


co 


0.8135 


1 = 1.229 - 1 = 0.229 


therefore 


m 


co 


™po 


= 4.366 


W CQ = 4.366(100) = 436.6 lb 
Wpo + W co = 100 + 436.6 = 536.6 lb 
Check by recalculating PcmSe* From equation (25), 


W 


co 


436.6 


Cm e N ttR 3 (55.5) (tt) (0.216) 

1N ma 7TK p 


= 11.58 lb/ft 3 


(Check is adequate: use 11.56 lb/ft 3 .) Check to see if bonding is necessary 

to prevent payload penetration 


R p\ 2 


m 


1 + 


co 

"Vo , 


= (0.286) 2 (5. 366) = 0.4386 


From figure 10 
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z s = 0.325 


Z max _ ^max _ ^Pmax _ 0.500 
z s Ps Ps 0 . 325 


1.538 > 1.000 


Therefore bonding is necessary to prevent payload penetration. 
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APPENDIX C 


SAMPLE CALCULATION FOR DESIGN BY SIMPLIFIED MODEL 
WITH PAYLOAD PENETRATION 


A sample calculation is given for an approximate design of a crushable 
casing under the assumptions that p c ^ = 0 , = 1 , and g = 0 (simplified 

model) together with Fpo( e ) = 1 for penetration. The design conditions are 
the following: 

U Q = 300 ft/sec 

ed = 0.7 

e m ~ 0 . 8 

n p ^ 2000 
Fmax 

a Q = 1200 psi for balsa- like material 

W D0 = m g = 450 lb 
P po e 

Rp = 1.2 ft 

Payload penetration permitted (unbonded) 

The last four conditions indicate that case 11 of table 1 is being calculated. 
From figure 5 for balsa- like material, 

SEA = 24,000 ft- lb/lb 


From equation ( 8 ), with o Q in psi and SEA in ft-lb/lb 


P g 
cm 6 e 


= £ m (144 ^ _ CO. 8 )C 144) Cl, 200) _ ... ,3 

SEA 24,000 • 0 ' 


From equation (22) with a Q in psi (for penetration) 


Pmax 


= ^ R p 2 Q44 g 0 ) ^ tt(i, 44 ) (144) Cl, 200) = 


m 


po & e 


450 


1,737 


For iteration .with figures 11 and 12, calculate 


1 / U ° 2 \ 1 C300) 2 

2e d yPmax®e R p J ~ X * 4 (1737) C32 . 17) ( 1 • 2 ) 


0.959 


49 


1 + 0.959 = 1.959 


(1 + 0 .959) 2 = 3.838 


SEA 


24,000 


n D e R 1737(0.8) (1.2) 
Fmax m p ^ 


= 14.385 


m, 


first 


co 


m 


=1.0 (guess) 


po 


first N 


= 1 + l '° 
mu 3.838 


= 0.5213 


From figure 11 for = 0.7 


From figure 12 


first = 0.8435 


r* P - 0.8435 _ n 

firSt R 1.959 0.4306 


first N mcr = 15.5 


m 1 _ 

second — — = _ — 1_ = 1.078 
m 14.385 
po 


second N 


mu 


1 + 1.078 
3.838 


= 0.5415 


From figure 11 for = 0.7 




second 

n ru = 

0.838 



second 

II 

0.838 _ 

1.959 " °' 4278 

From figure 

12 

second 

^ma 

15. 75 



third 

ii 

o o 

I 5 - 75 = 1 095 

14.385 ,U ^ 



third 

N 

mu 

1 + 1 . 095 n r . 

3.838 ’ °' 5460 

From figure 

11 for 

£ d = °- 7 
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From figure 12 

third N. 


ma 


Third N RU = 0.8368 


= 15.75 (same as second) 


R 


~~ = 0.4273 , R = 


1.2 


0.4273 


= 2.807 ft 


m 


co 


m 


= 1.095 


po 


W co = 450(1.095) = 493 lb 


Wp 0 + W co = 450 + 493 + 943 lb 


Check by recalculating p g 

cm e 


From equation (25) 


W 


p g = 

cm & e 


co 


493 




3 (15.75) (tt) (1. 728) “ 5.767 lb/ft 3 


(Check is adequate: use 5.76 lb/ft 3 ) 

Check that penetration occurs for this m /m and R r ,/R 

co po P 


r Rp\ 2 / m 

R - ) l 1 + m~~ ) = (0.4273) 2 (1 + 1.095) = 0.3825 
PO. 


From figure 10 


z = 0.257 
s 


From equation (5) with L = 0 


Z Pmax = e d[ 1 - R ~= °- 7 (°-5727) = 0.40089 


Zp max _ q Pma x _ 0.40089 
z s q s _ 0.257 


= 1.56 > 1.000 


Therefore penetration does occur 



APPENDIX D 

COMPUTER PROCEDURES FOR THE GOVERNING EQUATIONS 
WITH SPHERICAL GEOMETRY 


The computer procedures described for spherical geometry in this appendix 
form an essential part of the so-called "detailed" analytical model but can 
also be specialized for a variety of less detailed models. There are three 
computer procedures. 


BASIC COMPUTER PROCEDURE 


The first computer procedure described is the M basic M procedure. It 
cannot be regarded as a design since it contains no provision for automatic 
determination of crushable casing parameters to achieve a desired acceleration 
and/or a desired ratio of stroke to available stroke. The primary purpose of 
the basic procedure, then, is to check the adequacy of designs which have been 
determined by other means (for example, designs determined by the simplified 
model employing figs. 5 through 12). 

With symbols defined in the section on notation or in parentheses, the 
basic procedure contains the following steps: 

1. Input and print the following: case number, Wp Q , Rp, U Q , ej, cm, 

n md> g e * g L* g m ( value of g for equation of motion), g a (value of g for 
acceleration ratio) , (value of efc for equation of motion) , e]< a (value 

of efc for acceleration ratio), cks (value of for stress ratio), P c k m g e 

(value of Pc^ge ^ or equation of motion), P^gge (value of Pcfcge for 
acceleration ratio), P c ksge (value of P c k£e f° r stress ratio), Fp 0m (e) 

(value of F po (e) for equation of motion), F poa (e) (value of F po (e) for 
acceleration ratio) . Also input and print trial values of R and a 0 . 

Finally, input and print whether the payload is considered perfectly bonded 
or perfectly unbonded to the crushable material, whether the SEA is computed 
from o 0 or selected, whether the material is considered to be balsa- like or 
honeycomb-like (important only if SEA computed), and the value of SEA if 
selected. Note that PckinSe* PckaSe* °r Pcksge often is specified as p cm g e , 
to be calculated in step 2. Note also that Fp 0in (e) and Fp oa (e) is specified 

as a constant, usually 1.00, or as the integral of equation (A45) . 

2. Calculate and print the following constants: p p g e (from eq. (A43) 

with mp 0 = Wpo/ g e ) , SEA (from eq. (9) if calculated for balsa-like or 

eq. (10) if calculated for honeycomb- like) , P C mge (from eq. (8) with 

factor for dimensions), W co (from W co = (4/3)it(R 3 - R p 3 )p C mge)> W CQ + W po , 

Sm/^mdSe > ga/ n mdge> w o 2 (from the first of eqs . (A33)), K R (from the last of 
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eqs . (A33) with nip 0 = W po /g e and a factor for dimensions), and p pR g e (f rom 
eq. (A 35) with mpo = W po /g e ) . 

3. Generate and print values of F clin (from e q • (A36) with of 

step 1) and w 2 (from the equation of motion, eq. (A37) , with eq. (A34) and 

with cion, g m , and p c km of step 1) for selected values of z (z = q/R) . 
Terminate the integrations at the lowest positive value of w 2 when z is 
varying by increments of 0.0001. Print this value of z separately and label 
it z p max without penetration. Calculate L/R according to equation (5) with 

Z Pmax = ^Pmax^* P rint it separately, and label it L/R without penetration. 

4. Generate and print values of Fcls (from eq. (A36) with eks of 

step 1) and the stress ratio Op 0 k/ a 0 (from eq. (A38) with eq. (A34) and with 

and p c ks °f step 1) for selected values of z. With z varying by incre- 
ments of 0.0001, terminate the calculations at z = 1.00 or at the lowest 
positive value of 1 - ( a pok/ a o ) > whichever occurs at a lower z value. Print 
the termination value of z separately and label it z s . Calculate 

z D /z q . print it separately, and label it Zr, /z s without penetration. 

Fmax b r r ' Fmax * 

5. Generate and print values of F cla (from eq . (A36) for e^a of 

step 1) and deceleration ratio np/n mc i (from eq. (A32) with eq. (A34) and with 
e ka> ga> and Pcka from step 1) for selected values of z. Terminate the 
calculations at z p max without penetration. 

6. If the payload is bonded or if z n /z Q - 1.00 without penetration, 

nnax ^ 

terminate the program at step 5. If the payload is unbonded and if 
z Pma X / Zs > 1*00 without penetration, proceed into the penetration phase. For 

penetration, calculate the constant Kp according to the last of equations 
(A39) and calculate four initial conditions according to equations (A40) . 

7. Generate an d print v alues of y, e, dy/dx, and de/dx for selected 
values of x (x e t/n mc ig e /Rp) by integrating the two simultaneous ordinary 
differential equations (A41) , with auxiliary equations (A42) through (A45) , 
with initial conditions described in step 6, with F clm (y) determined as in 
step 3 when (Rp/R)y is substituted for z, with Fp 0m (e) being the integral 
of equation (A45) or a constant according to the specification in step 1, and 
with ekm> gm> and Pckm being the values of step 1. Terminate the integra- 
tion when dy/dx is zero to four or more decimal places and when successive 
values agree to four or more signifcant figures for the worst of y, e, and 
de/dx. The corresponding value of y is called y max , and the corresponding 
values of e and de/dx are initial conditions for the next phase of the 
problem. 

8. With the initial conditions just described, generate and print values 
of e and de/dx for selected values of x by integrating the first of equa- 
tions (A41) with d 2 y/dx 2 = 0. Terminate the integration when de/dx is 
zero to four or more decimal places and when successive values of e agree 

to four or more significant figures. The corresponding value of e is called 
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e max* Calculate Zp max from z pm&x = (R p /R) (y max + e kme m ax) > as deducible from 

equations (A39) with Zp = qp/R. Also calculate L/R from equation (5), and 

calculate Zr> / z q . Print z-n , L/R. and z-n /zc labeled "with 

, .. rmax s Rmax’ 5 Umax 5 

penetration. » 

9. Generate and print values of the deceleration ratio np/n mc j for 
selected values of x according to the first of equations (A41) with auxili- 
ary equations (A42) , (A43) , and (A45) . The quantity Fp 0a (e) is the integral 
of equation (A45) or a constant according to the specification in step 1; and 
e ka> &a* an d Pcka are th e values of step 1. Step 9 terminates the basic 
procedure when penetration is present. 


SEARCH FOR OVERALL RADIUS R 


The second computer procedure described is a design procedure in which 
o Q is assumed given but iterations are automatically performed to determine 
R for a desired ratio of stroke to available stroke. The desired ratio is 
unity in the present case, that is, L/R = 0; but the available stroke leaves 
a margin of safety based on the fictitious compacting strain e^. 

This procedure is based partially on a modified basic procedure. The 
modifications include the specification of nj es /n mc j (i.e., the desired value 
of n p max / n md) the determination and printout of z a (i.e., the value of 

z at n p = n pmax ) and 6 (i.e., ( n p max / n des) ' • 1116 quantity 3 is 

calculated as a measure of the acceleration discrepancy. 

The next step is to iterate the modified basic procedure just described 
in order to achieve a low value of L/R. As a start, the procedure is run for 
two values of R, the selection being based on the design charts, figures 6 
through 9, or any other analytical or experimental information suggesting low 
values of L/R. The results are labeled (L/R)! for R x and (L/R) 2 for R 2 . 

Then the iteration is based on successive straight lines of the form: 

L/R = aR + b (Dl) 

For the starting values (L/R)!, R x and (L/R) 2 , R 2 , the first pair of 
straight lines from equation (Dl) is 

(L/R)! = + b^ 

(L/R) 2 = a (1) R 2 + b O) 

and the computer determines a^ 1 -* and b*- 1 -* by solving the two simultaneous 
equations. From equation (Dl) for L/R = 0 (which is the desired value), the 

computer then determines R^ 1 ) as 
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I 


R ( 0 = -[b Cl Va (l) ] 


Then R^ is run through the modified basic procedure to determine (L/R)(l). 
This result is combined with the two starting values, and the two having the 
lowest absolute values of L/R are selected for a new pair of starting values. 
The process is repeated until 


|l/r| < 0.0005 

provided that each new value of | L/ R j is lower than at least one of its two 
starting values. 

If a new value of | L/R | is higher than either of its starting values, 
the computing machine is stopped, and new starting values must be selected. 
This did not happen for any of the cases reported herein, all of which con- 
verged rapidly (requiring an average of 0.17 min of execute time per case 
without penetration and 4.55 min for the detailed model with penetration). 

In contrast, there was one unrelated case that did not converge. For this 
case, however, convergence could not have been expected. It turned out that 
the given impact velocity was too high for a feasible energy absorbing design 
with the given material, and hence too high for a solution. 


SEARCH FOR OVERALL RADIUS R AND CRUSHING STRESS a 


The third computer procedure used herein is a search for R and a , which 
is simply an extension of the search for R just described (except that pene- 
tration is not included since eq. (28) determines a Q for penetration) . In 
the search for R and a Q , both R and a 0 are varied in an attempt not only to 
make L/R = 0 but also to make 6=0, where 6=0 when the maximum accelera- 
tion load factor n p max equals the desired load factor n^g. 

This additional requirement, 6=0, makes the iterations more complicated 
than before although the modified basic procedure to be iterated is the same. 
This time the starting values are found by running the basic procedure for 
three combinations of R and a 0 . Again, one of the starting values can be 
determined by the design charts for the simplified model, figures 6 through 9, 
together with modifications based on experience for other models and/or other 
preliminary information (such as analyses or experiments performed on similar 
configurations) . Experience indicates that the other two starting values 
should be small deviations from the first in which the higher values of o Q 
correspond to the lower values of R. The starting values are labeled (L/R)! 
and 6 i for R x and a Ql , (L/R ) 2 and 62 for R 2 ancl °c>2> anc * (L/R) 3 and 63 
for R 3 and cr o3 . The iterations are based on successive pairs of planar 
surfaces having the form 


L/R - aL 0 a o + 
6 = a 6a a o + 


&LR r + &L 

(D2) 

e 3rR + a B 

(D3) 
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For the three sets of starting values just labeled, the first group of 

three planes from equation (D2) is 

/ \ 




La 




. „(1). 


-I = a ia a o 2 + a TD R 2 + a 


,d), 

LR 


(l) 

L 


- a + a.™ 


a o3 + a D R 3 + a 


(i) 


^R/ 3 La ' “LR ’ “L 

and the computer solves the three simultaneous equations to determine a, 

and aRb Similarly, the first group of three planes from equation (D3) 
LK L 


CO 

La , 


is 


. a. (l) o„, ♦ a‘^1 


' a Bo'°ol * a"'Ri * a 


62 - a £ >002 ♦ aW « 2 ♦ a 6 


Cl) 

3 

(1) 


33 = a e ( a )a o3 + a g (l) R 3 + a 


CD 


and the computer determines af 1 ^, af ^ , and af 1 ^ . 

p0 pK p 

Then L/R and 8 are set equal to the desired value of 0 in equations 
(D2) and (D3) , and the coefficients just determined are substituted to give 

a Ci) a Ci) + a (i) R Ci) _ _ a Ci) 

a La °° a LR K a L 


a (1) R (1) - -a (l) 

3R R 3 

The computer then solves the two simultaneous equations and determines 


a (l) a (l) 
a 3a a ° 


(hopefully) an improved pair of parameters a*; 1 -* and R 


, 0 ) 


The parameters ao 1 '* and R*- 1 -* are introduced into the modified basic 

procedure to determine the corresponding values of (L/R)*- 1 -* and g^b This 
result is combined with the three sets of starting values, and the three hav- 
ing the lowest values of | L/R [ + | @ | are selected for new starting values. 
The process is repeated until 


|l/r| < 0.0005 
| 3 | < 0.0005 


with the requirement that each new value of | L/R | + |g| is lower than at 
least one of its three starting values. 
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The computing machine is stopped, analogously to the R search, if a 
new value of | L/R | + |$| is higher than any of its starting values; and new 
starting values must be selected [where the selection can often be facilitated 
by a plot having R and a Q as axes with values of | L/R | + |g| indicated by 
vectors or with values of L/R and 3 indicated by vectors) . In contrast to 
the R search, new starting values had to be selected fairly frequently for 
the search for R and cr 0 (specifically, for roughly one-fourth of the cases), 
despite the fact that convergence was rapid when it occurred (0.7 min of exe- 
cute time per case for the detailed model) . 


For all but one of the cases requiring a new set of starting values, 
convergence occurred for the second or third set. The exceptional case (not 
listed in table 1) was an attempt to do the detailed model by the R and o Q 
search so as to match an n p max °f 2000 selected for the simplified model. 

For the landing configuration under consideration, however, R searches for 


different values of o Q indicated that an 


as high as 2000 cannot be 


Pmax 

attained with the detailed model. Hence the iterative method would have been 
wrong if it had converged to an n p max of 2000 in a search for R and a 0 . 
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APPENDIX E 


EXACT INTEGRATION FOR A CLASS OF IMPACT PROBLEMS WITH 
VARIABLE MASS BUT WITHOUT PAYLOAD PENETRATION 


The class of impact problems to be considered is shown in figure 3(b). 

This shows a general vertically symmetrical landing geometry for zero shear 
resistance and uniform compacting strain with the impact to be described before 
payload penetration occurs (if any) . Equation (A19) defines the problem, and 
all assumptions leading to that equation are retained. 


One additional assumption is introduced, namely, a uniform and isotropic 
crushing stress (a = = constant) . Then equation (A19) becomes 


d(U 1 2 ) = g dq - ----- -■ 

po 


a k (A cih dq) 


m 


co 


- '‘■ckW 


CEl) 


where the equation has been multiplied by dq and where V c has been 
introduced according to the second of equations (A14) . 


The area Acih is the planar area of crushable material that is flush 
with the landing surface in figure 3(b), and V c is the volume of the dotted 
region in figure 3(b). Hence, 


Acih dq = dV c (E2) 

When equation (E2) is introduced into equation (El) , each term in the latter 
equation becomes a total differential, and the exact integration is 



gq 


+ 


£ k a k 


ck 


ic, g e 



+ m 


co 


“r)* 


(E3) 


where C is a constant of integration. This constant can be evaluated by 
noting that at q = 0, U = Uo, and V c = 0. When the evaluated C is intro- 
duced in equation (E3) , the result is 

1 2 112 ^ k k 

2 (u o - u ) = -gq + p — 

ck 


where V c has been replaced by 
equations (A14) . 


/' 


f 


-ck J A sh dh c /e k°v * m co> 
0 


•^sh dh c as ' t ^ ie second of 


(E4) 
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Equation (E4) gives U in terms of q and can be simplified if the gq 
term is neglected (as is justified for most impacts) . Then it is convenient 
to reintroduce V c and let it be the only variable on the right-hand side. 
Thus , 


1_ 

2 


(U 0 2 


- U2) 


¥k 

p ck 


log 


e 1 - [p . V /e, (m + m )] 
L ck c k^ po ccr J 


(E5) 


When U = 0, V c has reached its maximum value, V c . Thus, with the natural 
^ ^max 

log converted to an exponential: 


max 


e k (m p° * "co> 
p ck 


1 - e 


-( p ck u . z/! v k )' 


CE6) 


It is implicit in the derivation of equation (A13) that V c + V C1 = 
where Vci is the volume of m cl in figure 3(b); and the same relationship 
applies for the maximum volumes. Hence, if it is desired, equation (E6) can 
be written 


v c + V cl 
c max max 


V + m 


co 


ck 



ck U o /2e k CT k r 


(E7) 


Equation (E7) agrees exactly with equation (1-4) in appendix A of 
reference 2 when notational differences are accounted for. The development in 
reference 2 includes shear resistance implicitly (by showing uplift of the 
compacted region in a figure) but counteracts the inclusion of shear by assum- 
ing Aci to be a horizontal planar area and thereby rendering the shear 
stresses ineffective in energy absorption. Thus, with the results in agree- 
ment and the assumptions reconciled, the development of equations (E4) through 
(E7) becomes a partial check on the present methods as well as a set of exact 
results for a special case of variable mass in the absence of payload 
penetration . 
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TABLE 1.- NUMERICAL RESULTS FOR U 0 = 300 FT/SEC, e d = 0.7, e m = 0.8, AND np max < 2000 


Jq — JUU I l/OA-.Vj, C- d 

(a) Balsa-like and honeycomb -like material 



Rp» 

Ppge> 

lb/ft 3 


Material 

(a) 

Payload penetration 



Hp 

max 

V P si 

P 

H cnr 

! e . lb / ft3 ! 

Case 

ft 

Bonded 

Simplified 

Detailed 

Simplified 

Detailed 

Detailed 

Simplified! Detailed 

Detailed 

Simplified 

Detailed 

Detailed : 






00 

(c) 

(b) 

(c) 

Simplified 

(b) 

(c) 

Simplified 

tb) 

(c) 

Simplified 

1 

0.6 

110.5 

Yes 

H 

No 

No 

2000 

1999.7 

0.99985 

1080 

1 1395.0 

1.2917 

11.560 

13.294 

1.1500 

2 

.6 

110.5 

No 

H 

Yes 

No 

2000 

2000.4 

1.0002 

1228 

1395.8 

1.1366 

12.390 

13.299 

1.0734 

3 

.6 

110.5 i 

No 

B 

Yes 

Yes 

1954 

1954.8 

1.0004 

1200 

1200 

1 . 0000 

! 5.76 

5.7582 

.99969 

4 

.7 

69.6 

No 

H 

No 

No 

2000 

1999.9 

.99995 

1118 

1437.9 

1.2861 

11.82 

13.519 

1.1437 

5 

.7 

69.6 

No 

B 

Yes 

Yes 

2000 

1999.9 

.99995 

902.4 

1 902.4 

1.0000 

4.33 

4.3315 

1.0003 

6 

.8 

46.6 

No 

H 

No 

No 

2000 

1999.9 

.99995 

1170 

1491.4 

1.2747 

12.125 

| 13.796 

1.1378 

7 ' 

.9 

32.7 

No 

H 

No 

No 

1999.4 

2000.2 

1.0004 

1251.5 

( 1555.2 

1.2427 

12.518 

j 14.119 

1.1279 


SEA, ft-lb/1 

1 

b 

% /q s 

F max 5 

R, ft 



w c0 , lb 

case 

Simplified Detailed 

Detailed 

Simplified - Detailed 

Detailed 

Simplified 

Detailed 
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1 
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d Non- 
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exi stent 
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528.6 
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1.0017 1 

5 
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1 . 0000 
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2.3232 
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1.930 

1.9411 

1.0058 

125.0 
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1.0118 

6 
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1.1145 

dNon- 
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2.2595 

.95742 
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7 
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12,689 
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2.3883 
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lb 
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S36.6 

538.33 
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538.39 
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736.0 
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1.0058 

- . 000074 
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(a) H = honeycomb- like; B = balsa-like. 

(b) Design by simplified model employing figures 5-12 (except for cases 7 and 12 and the L/R column, which are machine computer results for 

simplified model) and assuming = 0, = 1, F^Ce) = 1» and g - 0. 


(c) Design by detailed model employing IBM 7094 solution of equations for P ck = P cm> = e^, F^ o (e) = 2 

(d) Over a hypothetical stroke as large as the outer radius, there is no force developed which is large enough to cause penetration. 
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NUMERICAL RESULTS FOR U Q = 300 FT/ SEC, ed = 0.7, e m = 0.8, AND n Pmax £ 2000 - Concluded 


(b) Balsa-like material with Wp 0 = 450 lb 
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No 
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No 

No 
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No 
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(c) 
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(b) 

(c) 
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8 
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3.456 

3.5288 

1.0211 

2.778 

2.9044 

1.0455 

329.0 

378.00 

1.1489 

9 

24,000 

24,007 

1.0003 

2.975 

2.7839 

.93573 

2.968 

2.9898 

1.0073 

598.2 

620.49 

1.0373 

10 

24,000 

24,000 

1.0000 

2.184 

2.1383 

.97908 

2.956 

3.0742 

1.0400 

391.0 

439.51 

1.1241 

11 

24,000 

24,007 

1.0003 

1.560 

1.3483 

. 86429 

2.807 

2.8457 

1,0138 

493.0 

514.17 

1.0429 

12 

24,000 

24,000 

1.0000 

1.4159 

1.2090 

.85387 

3.1545 
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1.0327 

460.76 
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1 . 11095 

13 

24,000 

24,000 

1.0000 

.730 

d Non- 
existent 


3.356 

3.4543 
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lb 
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Simplified 
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Detai led 
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Cc) 

Simplified 
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(c) 
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8 
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828.00 
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0.000886 

-.000162 
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9 
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1070. S 

1.0213 

-.001042 

.000172 

-.165 

10 

841.0 

889.51 

1.0577 

-.000525 

-.000063 

.12 

11 

943.0 

964.17 

1.0224 

-.000033 

.000290 

-8.79 

12 

910.76 

961.88 

1.0561 

.000329 

-.000057 

-.17 

13 

989.0 

1047.1 

1.0587 

-.000329 

.000383 
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(a) H = honeycomb-like; B = balsa-like. 

(b) Design by simplified model employing figures 5-12 (except for cases 7 and 12 and the L/R 

simplified model) and assuming = 0, e^ = 1, F^ o (e) = 1> and g = 0. 

(c) Design by detailed model employing IBM 7094 solution of equations for = p cm , = e m> 


colum^ which are machine computer results for 


F (e) 
P° 



(es p + l)s p dsp 
/2esp + e 2 + 1 


(d) Over a hypothetical stroke as large as the outet radius, there is no force developed which is large enough to 


cause penetration. 
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Stress 



(a) Stress-strain curves . 


Load 



Plate 

Crushed material 
Uncrushed crushable 

7777777777777777777777777777777777777777 ma ^ er,c, l 

(b) Crushing test with plate. 

Load 

Plunger 

Crushed material 

-Uncrushed crushable material 

7777777777777777777777777777777777777777 

(c) Crushing test with plunger. 

Figure 1.- Typical stress-strain curves and crushing tests. 
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Crushable 
material 
(mass=m C o) 


Landing 

777777777 / 

Surface 


m co m cr m pi 


i'/ 


(a) At the start of impact. (b) During impact, before 

payload penetration. 


(c) During impact, after 
payload penetration. 


Figure 2.- General vertically symmetrical landing geometry for zero shear resistance 





Payload 





qp-q 


7777777777777777777777777777777777Z75ffi7s 

Surface 


(a) At the start of impact. (b) During impact, before 

payload penetration. 


(c) During impact, after 
payload penetration. 


Figure 3.- General vertically symmetrical landing geometry for zero shear resistance and uniform 

compacting strain. 

















Hypothetical Nq 

r- Solid balsa 

cored balsa / 


SEA, ft-lb/lb 

0 

7.49xl0 6 /cr 0 ai 

□ 

(478.5) (cr o 0 - 4 ‘ 












Figure 11.- Mass and radius ratio design chart for simplified model with penetration when R p is given 



Figure 12.- Mass design chart for simplified model with penetration when Rp 
is given; also simple mass-volume relation. 


Crushable material weight, W c0 , 


Honeycomb-like material 
W po = IOOIb e m = 0.8 

n pmax " 2000 € d “ ® 

U 0 = 300 ft /sec 

Wco SEA 

o O Simplified model 1 ... . . .. 

□ a Detailed model j w,thout P enetratlon 



110.5 69.65 46.65 32.75 

Payload packaging density , p p g e , lb/ft 3 

(a) W co and SEA 

Figure 13.- Crushable casing properties and performance as functions of Rp 
and Ppg e for honeycomb- like material, a payload weight of 100 lb, and a 
payload maximum load factor of 2000. 
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Specific energy absorption , 
SEA, ftlb/lb 



Overall radius, 


See part (a) for design conditions 
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A 

Simplified model ] 
Detailed model J 

- without penetration 


A 

Simplified model 

with penetration 



Payload radius , R p , ft 


110.5 69.65 46.65 32.75 

Payload packaging density, yO p g e , Ib/ft 

(b) R and o Q 
Figure 13.- Continued. 
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Maximum crushing stress , 


Maximum value of payload displacement, q p q Pma ^ 

) 

q p at which penetration occurs or q s 

would occur if unbonded 


I 


See part (a) for design conditions 
Aim^e Qpmox /^s 

o Simplified model with penetration 



Payload radius , R p , f t 


110.5 69.65 46.65 32.75 

Payload packaging density, yO p g e , lb/ft 3 


Cc) ^Pmax /c ^s and Pcmge 
Figure 13.- Concluded. 
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Crushable material density, p crT] g e , Ib/ft 



Balsa-like material 
W p0 = 450 lb 

SEA = 24,000 ft lb /lb 
U 0 = 300 ft/sec 


m 


0.8 


cr Q = 800 psi 

Pc m g e ■ 3.84 lb/ft 3 

cr 0 = 1200 psi 

p g. =5.76 lb/ft 3 
r cm 3 e 


j- wi 
j- wi 


£ d =0.7 
n pmax * 2000 


without penetration 


with penetration 
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o 

x: 
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3 
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Undefined transitions for 
curves with unbonded payload 
(may require cr 0 < 1200 psi with 
penetration and cr Q > 800 psi 
without penetration) 



O Simplified model 
□ Detailed model 


} 


without penetration 


^ Simplified model \ with penetration 
i0) Detailed model J 

1 I I 1 

1.0 1.2 1.4 1.6 

Payload radius , R p , ft 

I I I I 


107.4 


62.15 


39.15 


26.23 


Payload packaging density, p g e , lb/ft' 


(a) W, 


co 


Figure 14.- Crushable casing properties and performance as functions of 
Ppge f° r balsa- like material and a payload weight of 450 lb. 


Lp and 
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Overall radius , 


See part (a) for design conditions 



107.4 62.15 39.15 26.23 

Payload packaging density, p g Q , Ib/ft 3 

00 R 

Figure 14.- Continued. 


81 








See part (a) for design conditions 

a„ Maximum value of payload displacement , q D 

“ Umnv » 


'Pmax 


q p at which penetration occurs or would occur if unbonded 



107.4 62.15 39.15 26.23 

3 

Payload packaging density, p p Q e , lb/ft 


(d ’ 'IPmax^s 
Figure 14.- Concluded. 
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